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ABSTRACT 


Most  gun  propellant  is  stored  at  depots  for  a  long  time 
before  it  is  used.  While  being  stored,  the  quality  of  the  gun 
propellant  may  deteriorate  and  become  unstable.  In  an  attempt 
to  avoid  disaster  due  to  use  of  unstable  gun  propellant, 
accurate  prediction  of  the  safe  shelf-life  of  gun  propellant 
is  necessary.  The  shelf-life  estimation  methods  used 
currently  for  a  group  of  similar  gun  propellant  lots  are 
based  on  a  fixed  effects  regression  model.  This  does  not  take 
into  consideration  the  fact  that  samples  from  the  same  lot  are 
more  similar  than  samples  between  lots.  To  capitalize  on  this 
lot-to-lot  variation  when  estimating  the  shelf-life,  first,  a 
random  effects  regression  model  is  developed.  Secondly,  a 
combined  mixed  effects  model  is  estimated.  The  estimated 
model  is  then  used  to  predict  not  only  the  shelf-life  of  a 
group  of  similar  lots  but  also  that  of  each  individual  lot  of 
5 " / 54  NACO  gun  propellant  stockpile.  The  results  indicate 
that,  first,  the  claimed  shelf -life  is  not  adequate  and 
requires  amendment.  Next,  the  group  shelf -life  estimated  can 
be  relatively  conservative  compared  to  the  individual  shelf- 
lives.  In  view  of  potential  opportunity  loss  due  to  safe 
individual  lots  being  discarded,  use  of  individual  shelf-life 
is  recommended. 
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EXECUTIVE  SUMMARY 


Most  gun  propellant  is  stored  at  depots  for  a  long  time 
before  it  is  used.  While  being  stored,  the  quality  of  the  gun 
propellant  may  deteriorate  and  become  unstable.  In  an  attempt 
to  avoid  disaster  due  to  use  of  unstable  gun  propellant, 
accurate  prediction  of  the  safe  shelf-life  of  gun  propellant 
is  necessary.  Before  1914,  real  time  storage  inspection  was 
used  to  estimate  the  safe  shelf-life  of  gun  propellant  lots 
produced  in  the  U.S.A.  In  an  effort  to  reduce  the  length  of 
time  required  for  real  time  storage  inspection,  an  Accelerated 
Aging  Test  (AAT)  was  introduced  at  Picatinny  Arsenal  in  1914. 
Under  AAT,  master  samples  taken  from  gun  propellant  lots  are 
sequentially  stored  in  heating  chambers  at  65.5  °C  to  measure 
the  red  nitrogen  oxide  fume  time.  Once  the  fume  time  is 
observed,  it  is  related  to  the  age  of  propellant  stockpiles. 

The  current  methods  to  estimate  the  shelf -life  of  a  group 
of  similar  gun  propellants  are  based  on  a  fixed  effects  linear 
regression  model  for  the  fume  time.  The  shortcomings  of  the 
current  estimation  methods  are  such  that  the  relationship 
between  the  fume  time  and  propellant  age  may  not  necessarily 
be  a  linear  one  and  the  fixed  effects  model  does  not  take  into 
account  the  potential  lot-to-lot  variation.  To  consider 
potential  lot-to-lot  variation,  a  random  effects  regression 
model  (RERM)  is  required. 

To  understand  the  state-of-art  estimation  of  the  safe 
shelf -life  of  gun  propellant,  four  technical  reports  are 
briefly  reviewed.  Furthermore,  a  preliminary  analysis  is 
conducted  to  understand  the  basic  characteristics  of  the 
master  samples  and  fleet  return  data  sets  of  5"/54  NACO  gun 
propellant  provided  by  Naval  Surface  Warfare  Center  (NSWC) . 
By  analyzing  the  fitted  models  and  validating  the  model 

vii 


assumptions,  the  standard  linear  regression  model  based  on  the 
reduced  data  set  is  selected  to  further  investigate  a  random 
effects  regression  model  analysis. 

To  consider  a  potential  lot-to-lot  variation  when 
estimating  the  shelf-life,  first  the  RERM  is  developed  using 
a  two-stage  analysis.  Next,  a  combined  mixed  effects  model  is 
estimated  by  employing  both  the  maximum  likelihood  (ML)  and 
restricted  maximum  likelihood  (REML)  methods.  These 
estimators  are  used  to  define  the  shelf-life  of  the  group  of 
similar  lots.  The  estimated  model  is  used  to  predict  not  only 
the  shelf-life  of  a  group  of  the  similar  lots  but  also  that 
for  each  individual  lot  by  a  shrinkage  procedure. 

The  results  indicate  that,  first,  the  claimed  shelf-life 
(35  years)  for  the  NACO  propellant  lots  is  not  adequate 
because  it  overestimates  the  shelf-life  for  NACO  propellant 
lots  by  almost  10  years.  Therefore,  it  requires  amendment. 
Next,  the  group  shelf -life  estimated  can  be  relatively 
conservative  compared  to  the  individual  shelf -lives.  In  view 
of  potential  opportunity  loss  due  to  individual  lots  discarded 
when  they  are  still  safe,  use  of  individual  shelf -life  is 
recommended  over  the  group  shelf-life  when  the  management  of 
individual  lots  is  expected  to  cost  less  than  that  of  bulk 
management . 

There  is  a  nonlinear  decreasing  trend  in  fume  time  in  some 
individual  lots  as  the  propellant  ages.  This  raises  questions 
about  the  approaches  that  use  traditional  linear  least  squares 
analysis  of  fume  time  data  to  predict  safe  shelf -life.  This 
thesis  recommends  that  a  nonlinear  regression  analysis 
approach  be  conducted  for  some  individual  propellant  lots. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Most  gun  ammunition  is  stored  at  depots  for  a  long  time 
before  being  used.  In  order  to  avoid  potential  disaster  due 
to  storing  unstable  gun  propellant,  the  shelf -life  of  a  gun 
propellant  stockpile  should  be  accurately  estimated. 
Typically,  the  shelf-life  of  gun  propellant  stockpile  is 
defined  as  the  age  at  which  5%  of  the  gun  propellant  is 
unstable.  The  estimated  shelf -life  can  also  be  used  to 
efficiently  minimize  the  ammunition  management  costs  for 
distribution,  maintenance,  consumption  and  disposal. 

In  the  early  1900' s,  real  time  storage  inspection  was 
used  to  estimate  the  safe  shelf-life  of  gun  propellant  lots 
produced  in  the  U.S.A.  However,  under  the  real  time  storage 
inspection,  it  took  a  long  time  to  find  visible  changes  in 
stability  of  propellent.  In  an  effort  to  reduce  the  length  of 
time  required  for  the  real  time  storage  inspection,  an 
Accelerated  Aging  Test  (AAT)  was  introduced  at  Picatinny 
Arsenal  in  1914  [Ref. 3].  For  AAT,  master  samples  are 
collected  for  testing  purposes  from  each  gun  propellant  lot 
produced  in  the  U.S.A.  over  the  past  sixty  years.  Each  master 
sample  consists  of  5  pounds  of  propellant  from  each  lot. 

Under  the  AAT,  forty- five  grams  of  propellant  are  taken 
from  each  master  sample  and  are  stored  in  heating  chambers  at 
65.5  °C  until  red  nitrogen  oxide  fumes  appear.  The  time  it 
takes  to  observe  such  oxide  fumes  is  called  the  fume  time. 
This  fume  time  is  recorded  with  the  corresponding  propellant 
age.  If  the  fume  time  falls  below  thirty  days,  the  propellant 
is  considered  unsafe  and  the  lot  is  recommended  to  be 
condemned.  A  fresh  batch  of  propellant  is  then  taken  from 
ambient  storage  and  the  AAT  is  repeated.  In  general,  it  is 
assumed  that  fume  time  decreases  as  propellants  age.  In  order 
to  validate  the  results  obtained  from  the  master  samples,  the 
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AAT  is  employed  on  the  fleet-return  data  of  gun  propellant 
lots  stored  in  the  active  depots. 

The  classical  methods  used  to  estimate  the  shelf-life  of 
a  group  of  similar  gun  propellants  is  based  on  a  fixed  effects 
linear  regression  model.  The  fume  time  data  are  aggregated 
over  a  group  of  similar  lots  and  a  linear  model  for  fume  time 
is  fitted  against  propellant  age.  The  safe  shelf -life  of 
current  methods  is  then  estimated  as  the  time  period  at  which 
a  95%  one  sided  lower  prediction  limit  for  the  fume  time  curve 
intersects  the  acceptable  lower  specification  level  of  thirty 
days  of  fume  time.  The  shortcomings  of  the  current  estimation 
method  are  such  that  the  relationship  between  the  fume  time 
and  propellant  age  may  not  necessarily  be  a  linear  one  and  the 
fixed  effects  model  does  not  take  into  account  the  potential 
lot-to-lot  variation.  Other  relationships  in  addition  to  the 
linear  model  should  be  examined.  In  order  to  consider 
potential  lot-to-lot  variations,  a  random  effects  model  is 
required. 

To  understand  the  state-of-art  related  to  the  estimation 
of  the  safe  shelf -life  of  gun  propellant,  four  technical 
reports  are  briefly  reviewed. 

B .  LITERATURE  REVIEW 

In  the  report  entitled  "Prediction  of  Safe  Life  of 
Propellants"  [Ref.  1],  the  nature  of  65.5°C  surveillance  test 
and  the  propellant  chemical  deterioration  are  briefly 
discussed.  In  this  study,  it  is  recognized  that  the 
measurement  of  residual  stabilizer  content  offers  the  best 
means  of  establishing  the  stability  potential  of  propellant. 
This  report  also  discusses  the  results  obtained  from  the 
standard  artillery  propellant  when  exposed  to  the  aging  test 
at  various  temperatures . 

The  measurements  of  residual  stabilizer  content  versus 
time  were  used  as  the  proxy  for  propellant  deterioration.  For 
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interpretation  of  this  deterioration  phenomena,  Berthelot's 
Law  was  employed.  It  demonstrated  that  a  family  of  straight 
lines  can  be  plotted  characterizing  the  length  of  time 
necessary  at  the  various  test  temperatures  to  obtain  a  given 
variation  of  stabilizer  content.  Finally,  by  establishing 
realistic  cut-off  points  regarding  stabilizer  content  in  a 
given  storage  temperature,  propellant  safe  life  was  estimated. 

Secondly,  the  report  entitled  "Type  Life  Program  History, 
Philosophy,  Accomplishment,  And  Future"  [Ref.  2]  introduces 
the  concept  of  the  type  life  test.  The  type  life  test  is 
derived  from  the  AAT  approach  and  takes  into  account  a  more 
realistic  temperature.  In  this  study,  accelerated  aging 
techniques  are  used  to  observe  the  aging  mechanisms  and 
characteristics  of  propulsion  ordnance  along  with  associated 
components.  They  are  also  use  to  determine  how  age  affects 
the  performance  of  the  units. 

The  type  life  temperature- time  profile  defined  in  this 
study  is  based  on  the  "hot-month"  concept,  and  differs  from 
the  AAT  and  its  use  of  constant  testing  temperature.  The 
concept  is  briefly  stated  as  follows : (1) propellants  are  a 
chemical  system; (2) chemical  systems  deteriorate  more  rapidly 
at  elevated  temperatures;  and  (3) in  a  normal  year's  time,  the 
deterioration  occurs  only  during  the  summer  months. 

Based  on  "hot-month"  concept,  the  engineers  establish  a 
temperature- time  profile  termed  the  "compressed-ambient  cycle" 
and  consisting  of  26  weeks.  This  profile  was  designed  to 
simulate  the  four  seasons  of  the  year  and  is  related  to  the 
aging  that  normally  occurs  during  magazine  storage  in  the 
U.S.A.  There  is  a  16-week  period  at  100  °F  (38 °C)  representing 
a  long  severe  summer  time,  two  3 -week  periods  at  70  °F 
representing  spring  and  autumn,  and  a  4 -week  period  at  40  °F 
representing  winter.  The  compressed- ambient  cycle  (26  weeks 
or  1/2  year)  simulates  one  year  of  magazine  storage;  therefore 
there  is  a  2:1  (i.e.52  weeks: 26  weeks)  aging  ratio.  The  main 
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purpose  of  this  report  was  to  evaluate  the  effect  of  the  loss 
of  nitroglycerin  from  the  propellant  into  the  inhibitor.  This 
loss  results  in  propellant  cracking  and/or  a  reduction  in 
propulsion  performance  due  to  a  loss  of  energy. 

Thirdly,  a  report  entitled  "Statistical  Analysis  of 
ARDEC'S  Fume  Data  Base"  [Ref.  3]  describes  the  procedure  for 
AAT  based  on  the  master  sample  surveillance  program  in  detail . 
Three  types  of  propellants  were  used  separately  based  on 
single  and  double  base  propellants.  The  single  and  double 
bases  are  terms  of  ammunition,  distinguished  by  the  methods  of 
ignition.  The  single  base  propellant  can  be  ignited  by  either 
mechanical  or  electrical  means.  The  double  base  propellant 
can  be  ignited  using  both  methods.  In  this  study,  fume  times 
are  plotted  against  sampling  times  for  each  propellant  form. 
Not  only  is  the  first  order  linear  relationship  examined  with 
this  approach  but  also  other  relationships  such  as  polynomial 
regressions,  logarithmic  and  exponential  transformation  are 
considered  to  establish  fume  time  behavior  over  the  life  of 
the  propellant. 

Finally,  a  study  entitled  "Long-Term  Stability  of  Navy 
Gun  Propellants"  [Ref.  4]  was  initiated  to  standardize  the 
requirements  of  safety  surveillance  testing  between  the  Army 
and  Navy.  This  was  necessary  because  current  Navy 
surveillance  testing  includes  a  65.5  °C  oven  fume  test  (AAT) 
on  single-base  propellants  and  double-base  propellants  with 
less  than  or  equal  to  10%  nitroglycerin  content. 
Additionally,  propellants  must  take  more  than  30  days  to  fume 
to  be  considered  stable.  In  contrast,  the  Army's  disposition 
criteria  is  based  on  a  stabilizer  analysis  of  less  than  0.2% 
nitroglycerin  content.  Because  a  substantial  quantity  of  Navy 
propellant  is  being  currently  stored  in  Army  depots,  non¬ 
standardization  complications  commonly  occur.  This  study  was 
tasked  by  Naval  Sea  Systems  Command (NSSC)  to  determine  the 
critical  level  of  propellant  stabilizer  for  selected 
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propellants  of  interest  to  the  Navy  and  to  the  Marine  Corps, 
to  develop  the  best  test  methods  to  evaluate  that  level  of 
criticality  for  Navy  propellants,  and  to  apply  the  data  for 
standardization. 

The  test  was  conducted  on  propellants  removed  from  loaded 
ammunition.  Each  sample  was  subjected  to  chemical,  kinetic, 
and  physical  properties  analyses.  As  a  result  of  these 
analyses,  it  was  recommended  to  continue  use  of  the  65.5  °C 
oven  fume  test  because  it  is  the  procedure  that  provides  the 
most  complete  coverage  of  propellant  in  storage  worldwide. 
The  Navy  has  a  large  historical  database  of  fume  data,  and 
statistical  analysis  of  the  data  provides  an  inexpensive, 
comprehensive  means  of  insuring  that  Navy- developed  propellant 
is  safe  for  worldwide  storage. 

In  summary,  the  master  sample  surveillance  test  is 
recommended  for  continuous  use  to  determine  safety  of  the  gun 
propellant  in  storage.  It  is  proposed  that  the  relationship 
between  the  fume  time  and  propellant  age  may  not  be  limited  to 
a  linear  relationship.  The  management  of  data  structure  and 
the  data  transformation  (e.g.  logarithmic  or  reciprocal)  is 
required.  Additionally,  no  potential  lot-to-lot  variation  in 
shelf-life  estimation  was  examined  in  the  four  technical 
reports  summarized  here. 

C.  THE  SCOPE  OF  THESIS 

In  this  thesis,  the  master  sample  surveillance  test  data 
set  is  first  examined  using  several  fixed  effects  regression 
models  in  an  attempt  to  find  the  relationship  between  the  fume 
time  and  propellant  age.  Then,  a  random  effects  regression 
model  is  introduced  to  take  into  account  a  lot-to-lot 
variation  in  fitting  the  relationship  between  the  fume  time 
and  propellant  age.  Finally,  results  are  discussed  and 
recommendations  are  made. 
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II.  FIXED  EFFECTS  REGRESSION  MODEL 


A.  DATA  DESCRIPTION 

A  preliminary  analysis  is  conducted  to  understand  the 
basic  characteristics  of  the  master  sample  data  set  provided 
by  Naval  Surface  Warfare  Center  (NSWC)  .  The  data  set  was 
collected  based  on  the  oven  test  data  consisting  of  28  lots  of 
master  samples  of  5"/54  NACO  gun  propellant.  They  were 
manufactured  at  Badger  Army  Ammunition  Plant  (BAAP) ,  Baraboo, 
Wisconsin  between  1972  and  1974  .  The  NACO  is  a  generic  term 
of  NAvy  COol  representing  a  family  of  cool -burning,  single¬ 
base  propellants.  These  propellants  are  made  from  low- 
nitration  nitrocellulose  with  a  coolant,  a  stabilizer,  a 
decoppering  agent,  and  a  flash  suppressant  added. 

First  by  subtracting  the  manufacturing  date  (see  Appendix 
A)  from  the  date  when  a  45g  propellant  batch  was  sent  to  the 
heating  chamber,  the  age  is  computed  for  each  propellant  lot. 
The  data  of  the  fume  time  and  propellant  age  is  shown  in 
computer  program  list  in  Appendix  B.  Both  descriptive 
statistics  of  fume  time  and  propellant  age  are  displayed  in 
Table  1.  In  summary,  221  observations  of  fume  time  and 
propellant  age  had  been  taken  from  28  propellant  lots.  The 
maximum  (1492  days)  and  minimum  (19  days)  of  fume  time  are 
observed  during  6813  days  of  propellant  age.  Next,  in  an 
attempt  to  find  the  relationship  between  the  fume  time  and 
propellant  age,  the  fume  time  is  plotted  against  the 
propellant  age  in  Figure  C.l.  It  is  observed  that  there  is  an 
unexpected  increasing  trend  in  fume  time  over  propellant  age 
up  to  1100  days.  This  unexpected  pattern  may  be  due  to  the 
type  of  ammunition,  shape  of  propellant  container  or  caused  by 
some  unknown  chemical  reaction.  [Ref.  3] 
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Descriptive  Statistics  For  Raw  Data 

manufacture  date 

min  1972  max  1974 

number  of  lot 

28 

number  of  repeated 

measurements 

221 

fume  time 

in  days 

mean 

802.61 

std 

222 . 159 

max 

1492 

min 

19 

propellant 
age  in  days 

mean 

3488 . 824 

std 

2020 . 824 

max 

6813 

min 

0 

Table  1.  Descriptive  Statistics  For  Raw  Data 


B.  STANDARD  LINEAR  REGRESSION  MODEL  I 

Regardless  of  such  unexpected  pattern,  according  to  the 
current  practice,  standard  linear  regression  model  I,  is 
applied  to  fit  the  fume  time  (Y)  against  propellant  age  (X) : 

Y  =  Bq  +  Bx  X  +  e  (1) 


where  B0  and  Bj  represent  the  expected  initial  condition  and 
deterioration  rate  of  stability  of  gun  propellant,  and  e  is 
the  random  error  which  follows  a  Normal  (0,t2)  distribution. 
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The  predicted  value  ( ? J  at  X=X»  is 

Y*  =  4  +  4  Z*  (2) 

where  B0  and  Bx  are  the  ordinary  least  square  estimators  of  B0 

and  Bj .  Subsequently  a  90%  prediction  interval  for  a  fume 
time  at  given  propellant  age  (X*)  is 

Yt  ±  t  {  0.95  ,  N- 2  )  SEPRED  (  Y,  \  X,  )  ,  (3) 

where  N  is  the  total  number  of  observations,  t(0.95,N-2)  is 
the  95th  percentile  of  the  t  distribution  with  N-2  degrees  of 
freedom.  SEPRED  (Yt  \Xt)  is  the  standard  error  of  predicted  Yt 
at  X=X* : 

SEPRED  (  Yt  |  X,  )  =  T  [1  +  4  +  _}  — - —  ]  7  .  (4) 

N  £  (  "  X  )  2 

Based  on  this  information,  the  group  shelf-life  (tsl)  can  be 
estimated  by  checking  the  time  period  when  the  95%  one-sided 
lower  prediction  limit  intersects  30  days  of  fume  time: 

30  =  B0  +  4  tsJ  -  t  (  0.95  ,  N-2  )  SEPRED  (  Y,  \  tsl  )  .  (5) 

Next,  actual  data  obtained  from  the  NACO  propellant 
surveillance  test  is  used  to  fit  the  model  and  to  estimate  the 
shelf -life.  The  fitted  model  based  on  (2)  is  given  in  Table 
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Standard  Linear  Regression  Model  I 

/V 

B0 

1045.510 

/V 

B1 

-0 . 070 

SE(B0) 

23 . 171 

SEiBj 

0.006 

R-square 

0.4009 

lack  of  fit  test 

degrees  of  freedom:  200,19; 

P  value  =  0 

Normal  distribution  test 

P  value  =  0 . 004 

(Kolmogorov-Smirnov  Test) 

not  Normal 

estimated  shelf-life 

10400  days  (28.49  years) 

Table  2.  Fitted  Standard  Linear  Regression  Model  I 


2  and  displayed  as  a  solid  line  in  Figure  C.l.  A  90%  two- 
sided  prediction  interval  (solid  curves)  is  overlaid  in  Figure 
C.l.  The  resultant  shelf -life  based  on  (5)  turns  out  to  be 
10400  days  (28.49  years)  as  shown  in  Figure  C.2. 

The  residual  analysis  [ri  =  Y1  -  Y±  ,  for  i=  1 ...  N)  for 

standard  linear  regression  model  I  is  shown  in  Figure  C.3. 
Residual  pattern  indicates  nonlinearity  and  nonconstant 
variance.  Also,  the  Kolmogorov- Smirnov  test  statistic  shown 
in  Table  2  indicates  that  the  residuals  do  not  follow  a  normal 
distribution . 

For  further  data  analysis,  the  data  set  that  covers  from 
zero  to  1100  days  of  propellant  age  is  deleted,  because  the 
pattern  is  unstable  during  this  early  stage  and  this  pattern 
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is  not  expected  to  last  longer  than  1100  days.  Furthermore, 
after  1100  days  of  propellant  age,  the  fume  time  tends  to 
decrease  as  propellants  age  and  this  pattern  is  more 
consistent.  Therefore,  the  use  of  the  remaining  data  set 
appears  to  be  reasonable  for  the  shelf -life  estimation.  The 
following  three  regression  models  are  employed  (1)  standard 
linear  regression  model  II  (delete  observations  that 
correspond  to  propellant  age  below  1100  days  and  use  standard 
linear  regression  model  I);  (2)  linearizing  transformation 
model  I  (transform  propellant  age  to  1/  propellant  age) ;  and 
(3)  linearizing  transformation  model  II  (transform  fume  time 
to  log (fume  time)).  Graphical  method  precedes  statistical 
analysis . 

C.  STANDARD  LINEAR  REGRESSION  MODEL  II 

The  standard  linear  regression  model  II  is  essentially 
the  same  as  the  standard  linear  regression  model  1(1)  . 
However,  parameters  are  estimated  based  on  a  reduced  data  set. 
After  deleting  the  early  age,  the  fume  time  is  plotted 
against  propellant  age  in  Figure  C.4.  The  results  of  model 
fit  are  shown  in  Table  3 .  The  predictive  value  (middle  solid 
line)  and  a  90%  two-sided  prediction  interval  for  a  fume  time 
(solid  curves)  are  overlaid  in  Figure  C.4. 

The  residual  analysis  for  the  standard  linear  regression 
model  II  is  shown  in  Figure  C.5.  In  general,  the  plot  of 
residuals  against  fitted  values  over  fume  time  between  zero 
day  and  750  days  (propellant  age  above  4500  days)  does  not 
show  any  systematic  features  and  they  appear  to  have  a  common 
variance.  But,  a  minor  nonrandom  pattern  is  observed  when  the 
fume  time  exceeds  750  days  (propellant  age  below  4500  days) . 
Also,  based  on  the  p- value  associated  with  Kolmogorov- Smirnov 
test  shown  in  Table  3,  we  can  conclude  the  residual  of  this 
model  follows  a  normal  distribution. 


11 


Standard  Linear  Regression  Model  II 

Bq 

1388.452 

/\ 

Bi 

-0 . 141 

SE(B0) 

20 . 612 

SE  ( Bx ) 

0 . 005 

R- square 

0.838 

lack  of  fit  test 

degrees  of  freedom:  167,11; 

P  value  =  0.001 

Normal  distribution  test 

P  value  =  0.458 

(Kolmogorov- Smirnov  Test) 

Normal 

estimated  shelf-life 

8469  days  (23.20  years) 

Table  3.  Fitted  Standard  Linear  Regression  Model  II 


By  checking  the  level  of  95%  one-sided  lower  prediction 
limit  for  a  fume  time  that  intersects  3  0  days  of  fume  time  in 
Figure  C.6,  the  estimated  group  shelf -life  turns  out  to  be 
8469  days  (23.20  years). 

D.  LINEARIZING  TRANSFORMATION  MODEL  I 

In  linearizing  transformation  model  I,  it  is  assumed  that 
the  fume  time  is  linearly  related  to  reciprocal  of  propellant 
age : 

r  =  +  1  )  +  e  (6) 

The  fume  time  is  plotted  against  propellant  age  in  Figure  C.7. 
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The  results  of  a  model  fit  are  shown  in  Table  4. 


The 


predicted  value  for  the  fume  time 


y.  =  Bn 


/\  , 

B1  ( 


X, 


is  shown 


Linearizing  Transformation  model  I 

/V 

B0 

411.738 

1315827.380 

SE(B0) 

17.420 

SE  { Bx ) 

51308.706 

R- square 

0 . 787 

lack  of  fit  test 

degrees  of  freedom:  167,11; 

P  value  =  0 

Normal  distribution  test 

P  value  =  0 . 02948 

|  (Kolmogorov-Smirnov  Test)  i 

not  Normal 

estimated  shelf-life 

00 

Table  4 .  Fitted  Linearizing  Transformation  Model  I 

as  a  middle  curve  and  a  90%  two-sided  prediction  interval  for 
a  fume  time  are  overlaid  in  Figure  C.7. 

It  is  interesting  to  note  that  the  95%  one-sided  lower 
prediction  limit  for  predictive  fume  time  converges  to  around 
400  days  as  the  propellant  age  continuously  increases. 
Therefore,  the  safe  shelf-life  of  gun  propellant  based  on 
linearizing  transformation  model  I  cannot  be  estimated  or  can 
be  considered  as  infinity. 

The  residual  analysis  for  the  linearizing  transformation 
model  I  is  shown  in  Figure  C.8.  In  general,  the  plot  of 
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residuals  against  fitted  values  over  fume  time  between  400 
days  and  800  days  (propellant  age  above  4000  days)  does  not 
show  any  systematic  features  and  they  appear  to  have  a  common 
variance.  But  a  nonrandom  pattern  is  observed  when  the  fume 
time  exceeds  800  days  (propellant  age  below  4000  days) .  This 
phenomenon  is  similar  to  the  standard  linear  regression  model 
II  (see  Figure  C.5).  Also,  based  on  the  p-value  associated 
with  Kolmogorov -Smirnov  test  shown  in  Table  4,  we  can  conclude 
that  the  residual  of  this  model  does  not  follow  a  normal 
distribution . 

E.  LINEARIZING  TRANSFORMATION  MODEL  II 

In  linearizing  transformation  model  II,  it  is  assumed 
that  the  propellant  age  is  linearly  related  to  log  fume  time: 

LOG  ( K)  -  B0  +  Bx  X  +  e.  (7) 


The  fume  time  is  plotted  against  propellant  age  in  Figure  C.9. 
The  results  of  a  model  fit  are  shown  in  Table  5.  The 
predicted  value  for  the  fume  time  Y„  -  EXP  ( B0  +  B xX,)  and  a  90% 
two-sided  prediction  interval  for  a  fume  time  are  overlaid  in 
Figure  C.9. 

The  residual  analysis  for  the  linearizing  transformation 
model  II  is  shown  in  Figure  C.10.  In  general,  the  plot  of 
residuals  against  fitted  values  over  fume  time  between  zero 
day  and  6.6  days  (propellant  age  above  4  0  00  days)  does  not 
show  any  systematic  features  and  they  appear  to  have  a  common 
variance.  But  a  nonrandom  pattern  is  observed  when  the  fume 
time  exceeds  6.6  days  (propellant  age  below  4000  days) .  This 
phenomenon  is  also  similar  to  the  previous  two  models  (see 
Figure  C.5  and  Figure  C.8).  Also,  based  on  the  p-value 
associated  with  Kolmogorov-Smirnov  test  shown  in  Table  5,  we 
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Linearizing  Transformation  model  II 

B0 

7.38 

-1.79  E-4 

SE(B0) 

5.75  E-2 

SE  ( Bx ) 

1.29  E-5 

R- square 

0.520 

lack  of  fit  test 

degrees  of  freedom:  167,11; 

P  value  =  0 

Normal  distribution  test 

P  value  =  1 . 0E-11 

(Kolmogorov -Smirnov  Test) 

not  Normal 

estimated  shelf-life 

17400  days  (47.67  years) 

Table  5.  Fitted  Linearizing  Transformation  Model  II 


can  conclude  that  the  residual  of  this  model  does  not  follow 
a  normal  distribution. 

By  checking  the  level  of  the  95%  one-sided  lower 
prediction  limit  of  a  fume  time  that  intersects  30  days  of 
fume  time  in  Figure  C.ll,  the  estimated  group  shelf -life  turns 
out  to  be  17400  days  (47.67  years) . 

In  summary,  by  analyzing  the  fitted  models  and 
validating  the  model  assumptions  of  the  various  fixed  effects 
regression  models,  the  effective  distribution  data  model 
appears  to  be  the  standard  linear  regression  model  II  based  on 
the  reduced  data.  This  model  is  compared  to  the  model  fit 
based  on  the  fleet  return  data. 
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F.  STANDARD  REGRESSION  MODEL  BASED  ON  FLEET  RETURN  DATA 

This  model  is  based  on  the  oven  test  data  consisting  of 
28  lots  of  fleet  return  data  of  5"/54  NACO  gun  propellant. 
The  descriptive  statistics  of  fume  time  and  propellant  age  are 
displayed  in  Table  6.  In  summary,  204  observations  of  fume 


Descriptive  Statistics 

manufacture  date 

min  1972  max  1974 

number  of  lot 

28 

number  of  repeated 

measurements 

202 

fume  time 

in  days 

mean 

702.12 

std 

85.46 

max 

993 

min 

479 

propellant 
age  in  days 

mean 

2822 

std 

1388.8 

max 

5074 . 5 

min 

361 

Table  6 .  Descriptive  Statistics  For  Fleet  Return  Data 


time  and  propellant  age  had  been  taken  from  28  propellant 
lots.  The  maximum  (993  days)  and  minimum  (479  days)  of  fume 
time  are  observed  during  5074.5  days  of  propellant  age.  The 
fume  time  is  plotted  against  propellant  age  in  Figure  C.12. 
The  results  of  model  fit  which  is  based  on  standard  linear 
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Standard  Linear  I 

Degression  Model 

/V 

823.375 

-0 . 042 

SE(B0) 

9 . 793 

SEiBj 

0 . 003 

R- square 

0.422 

lack  of  fit  test 

degrees  of  freedom:  127,73; 

P  value  =  0.0000 

Normal  distribution  test 

(Kolmogorov- Smirnov) 

P  value  =  0.0040 

not  Normal 

estimated  shelf-life 

15700  days  (43  years) 

Table  7.  Fitted  standard  Linear  Regression  Model 
Based  on  Fleet  Return  Data 

regression  model  I  are  shown  in  Table  7.  The  predicted  value 
(middle  solid  line)  and  a  90%  prediction  interval  for  a  fume 
time  (solid  curves)  are  overlaid  in  Figure  C.12.  The 
resultant  shelf-life  based  on  (5)  of  this  model  turns  out  to 
be  15700  days  (43  years)  as  shown  in  Figure  C.13. 

The  residual  analysis  for  this  model  is  shown  in  Figure 
C . 14 .  There  is  a  right -opening  megaphone  pattern  in  Figure 
C . 14 .  This  residual  pattern  observed  shows  nonconstant 
variance.  Also,  the  Kolmogorov -Smirnov  test  statistic  shown 
in  Table  7,  indicates  that  the  residuals  of  this  model  do  not 
follow  a  normal  distribution.  Two  remedies  for  nonconstant 
variance  (Weisberg  [Ref.  8])  are  made  using  weighted  least 
square  estimation  and  power  transformation: 
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1.  Weighted  Least  Square  Model 

The  first  remedy  is  the  weighted  least  square  model.  The 
variance  of  the  random  error  in  the  standard  linear  regression 
model  I  (1)  appears  to  increase  as  X  increases,  that  is,  r=KX 
where  K  is  a  positive  constant.  Therefore,  the  model  (1)  is 
modified  as  follows: 


—  =  —  +  B.  +  — 
X  X  1  X' 


(8) 


The  residual  analysis  for  this  model  is  shown  in  Figure  C.15. 
It  still  shows  a  right -opening  megaphone  pattern  in  Figure 
C.15  and  the  suggested  weighted  least  square  method  is  not 
adequate . 

2.  Power  Transformation  Model 

The  second  remedy  is  to  transform  the  response  Y  via  a 
variance  stabilizing  transformation.  By  checking  the  residual 
pattern  of  several  transformation  models,  the  following 
transformation  model  appears  to  satisfy  the  assumption  of 
common  variance  (see  Figure  C.16) : 

Y~3  =  B0  +  B1  X  +  e.  (9) 


The  fume  time  is  plotted  against  propellant  age  in  Figure 
C.17.  The  results  of  a  model  fit  are  shown  in  Table  8.  The 

_JL 

prediction  fume  time  Y  =  (B0  +  Bx  Xt)  3  and  a  90%  prediction 

interval  for  a  fume  time  are  overlaid  in  Figure  C.17. 

The  p-value  associated  with  Kolmogorov- Smirnov  test  of 
residual  of  model  (9)  indicates  that  the  residual  of  this 
model  follows  a  normal  distribution. 
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Power  Transformation  model 

4 

1.75  E-9 

A 

4.87  E- 13 

SE(B0) 

1.26  E- 10 

SEiBj 

4  E-14 

R- square 

0.425 

lack  of  fit  test 

degrees  of  freedom:  127,73; 

P  value  =  0.001 

Normal  distribution  test 

P  value  =  0.18964 

( Kolmogorov- Smirnov) 

Normal 

estimated  shelf-life 

00 

Table  8 .  Fitted  Power  Transformation  Model 

The  shelf-life  is  estimated  based  on  the  following 
formula : 


-1  (10) 

30  =  [B0  +  B±  tsl  -  t(  0.95  ,  N-2  )  SEPRED (  Yt  \  tsl  )]  3  . 


The  estimated  group  shelf-life  turns  out  to  be  infinite. 

This  shelf-life  estimation  based  on  fixed  effects 
regression  models  does  not  take  into  account  the  lot-to-lot 
variation.  Apparently  potential  lot-to-lot  variation  is 
supported  by  individually  fitted  models  based  on  the  Standard 
Linear  Regression  Model  I,  II  and  fleet  return  data  as  shown 
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in  Figures  C.18,  C.19,  C.20.  Additionally,  this  variation  can 
be  seen  by  checking  linearizing  transformation  model  I  after 
smoothing  the  data  following  Lowess  procedure  [Ref.  9]  as 
shown  in  Figure  C.21. 

In  the  next  chapter,  standard  linear  regression  model  II 
will  be  further  investigated  for  a  random  effects  regression 
model  analysis. 
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III.  RANDOM  EFFECTS  REGRESSION  MODEL 


A.  INTRODUCTION 

The  current  methods  used  to  estimate  the  shelf-life  of 
gun  propellant  is  based  on  the  standard  linear  regression 
model  I.  In  this  fixed  effects  linear  regression  model,  the 
expected  initial  fume  time  (B0)  and  deterioration  rate  (Bj)  are 
assumed  to  be  constant  for  all  gun  propellant  lots.  However, 
in  practice,  the  performance  varies  among  lots  even  though  a 
group  of  lots  is  manufactured  based  on  the  same  technique. 
This  kind  of  phenomenon  can  be  observed  in  the  individually 
fitted  models  in  Figure  C.18  to  C.21.  The  lot-to-lot 
variation  is  observed  not  only  in  the  intercepts  (expected 
initial  condition)  but  also  in  slopes  (deterioration  rate) 
among  a  group  of  similar  lots.  The  current  methods  does  not 
take  into  account  this  kind  of  constant.  In  order  to  consider 
potential  lot-to-lot  variation  within  the  group,  a  random 
effects  (RERM)  (or  mixed  effects)  regression  model  is 
necessary. 

In  this  chapter,  first,  the  RERM  will  be  introduced  based 
on  a  two-stage  analysis.  Next,  a  combined  two-stage  model 
(mixed  effects  model)  will  be  fitted  based  on  the  reduced  data 
of  the  master  samples  provided  by  NSWC .  In  order  to  obtain 
the  estimators  of  unknown  parameters  in  mixed  effects  model, 
the  maximum  likelihood  estimation  (MLE)  and  restricted  maximum 
likelihood  (REML)  will  be  employed.  These  estimators  are, 
then,  used  to  defined  the  shelf-life  of  the  group  of  similar 
lots  that  takes  into  account  the  lot-to-lot  variation. 
Furthermore,  a  shrinkage  procedure  will  also  be  introduced  to 
estimate  the  shelf -life  of  an  individual  lot. 
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B .  MODEL 

1.  Within-Lot  Model 

In  order  to  estimate  the  deteriorating  patterns  of  gun 
propellant  lots  over  storage  time,  each  45g  batch  is  taken 
from  5  pound  master  sample  of  lot  i  and  its  ageing  is 
accelerated  at  the  heating  chamber  at  65.5  °C.  The  fume  time 
(yy)  of  lot  i  (i=l...N)  is  recorded  with  the  corresponding 
propellant  age  (ty)  (j=l...n;)  whenever  the  red  nitrogen  oxide 
fume  appears.  Once  repeated  measurements  are  obtained,  we 
attempt  to  find  the  relationship  between  the  fume  time  and 
propellant  age.  The  following  within-lot  model  is  used  to 
describe  the  relationship  between  the  fume  time  and  the  age  of 
a  lot : 

For  i  =  l...N  and  j=l...n; 


y. 
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tij  +  €ii 
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where  aoj  and  a,,  represents  the  expected  initial  condition  and 
deterioration  rate  of  the  stability  of  propellant  lot  i,  and 
£  y  is  the  random  error  which  is  assumed  to  follow  a 
independent  N(0,  r2)  . 

2 .  Between-Lot  Model 

In  order  to  take  into  account  the  lot-to-lot  variation 
among  intercepts  (expected  initial  condition)  and  slopes 
(deterioration  rate)  of  a  group  of  lots,  we  assume  the 
following  between-lot  model: 


Oi 

=  w± 

Yo 

+  Soi 

1  i 

=  wi 

Yi 

+  »ii 

(12) 
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where  Wj  is  a  1  x  k  vector  representing  values  of  K  factors 
that  may  cause  varying  a0j  and  au  .  The  y0  and  yl  are  k  x  1 
vectors  of  regression  coefficients  for  aoi  and  c^, 
respectively.  The  <50i  and  5n  are  assumed  to  follow  independent 


„  ol  Oqi 

N(0,£),  where  L  =  .  Note  that  for  AAT,  the  master 

oh  0i 


samples  are  put  in  the  constant  temperature  condition  (65.5 
°C)  and  W;  in  this  case  is  1. 

When  a0i  and  an  in  (11)  are  substituted  with  those  in 
(12)  ,  a  combined  mixed  effects  model  and  is  formulated  as 
follows : 

For  i  =  l...N  and  j=l...nj, 


This  mixed  effects  model  can  be  re-written  using  matrix 
notation : 


y  =  Xp  +  Zv  +  e  ,  (14) 


where  y  is  a  .n^xl  vector  of  y^;  /3  is  a  2Kxl  coefficient  vector(y0  ,  Yi)7 

N 

associated  with  a  V}  nix2Jf  matrix  of  X  consisting  of 
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(Wj,  ,  Wit^)  ;  v  is  a  2n.jXl  vector  of  (80i  ,  b1±y  associated  with 

N 

a  known  V  nix.2ni  block  diagonal  matrix  Z  of  (1  ,  t1;j)  ,  and  e 

fz i 

N 

is  an  £  iijXl  vector  of  random  error.  Assume  that  v  and  e 

are  uncorrelated  and  have  expections  0  and  variances  G  and  R, 
respectively,  where  both  G  and  R  are  nonsingular.  applying 
these  assumptions  to  (14)  leads  to 


E(y)  =  xp 

Var{y)  =  V  =  ZGZ1  +  R  , 


(15) 


where  G  is  a  2njX2ni  block  diagonal  matrix  consisting  of 
and  R  =  Var(e)  =  t2I» 

When  G  and  R  (and  hence  V)  are  known,  estimates  of  both 
(3  and  the  realized  value  of  v  are 


{X'V^X)  '1X/V1y 


v=  [vx,  .  .  .  vN]  '=  [$01, 6] 


respectively,  where  P  is  the  best  linear  unbiased  estimator 
(BLUE)  of  13,  and  v  is  the  best  linear  unbiased  predictor 
(BLUP)  of  v.  They  can  also  be  obtained  from  Henderson's 
(1984)  mixed  model  equation  [Ref.  10]  : 
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'  X'R-'X  X'R~xZ  1  [  p 

Z’R^X  Z'R^Z+G'1  J  [  $ 

where  p  =  [?0  ,  fj'  and  $  =  [801  ,  6^,  .  .  ,§0N  ,  81A]'.  When  G  and  R 

(and  hence  V)  are  unknown,  the  estimator  of  V  can  be  obtained 
by  maximizing  the  logarithm  of  likelihood  function  [Ref.  10] 


log  L=- -i  N  log27t--ilog|v|--i  (y-Xfi) 1  V  1  (y-x P)  .  (18) 

The  ML  estimators  V  and  P  can  be  obtained  by  solving  the 
following  equations  simultaneously 

X'  V^1  X  p  =  X1  tf"1  y  <19) 

and 

tr  ( V"1ZiZ/1)  =  (y-xp)  '1F1ziZ,iQ~1  (y-Xp)  for  i=0 , 1 ,  .  .  . N,  (20) 

N 

where  Z{  is  a  S' nix2ni  block  diagonal  matrix  of  ( 1 ,  t y )  and 
T=i 

tr  ( V'"1ZiZJ/)  is  trace  of  the  ^  nix^2  ni  matrix  {V~1ZiZl1)  that 

2=1  i=l 

denotes  the  sum  of  the  diagonal  entries.  An  algebraically 
simpler  expression  for  (20)  is  derived  by  defining 
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p  =  V~ 1  -  V-1X{XIV~1X)~1XIV-1 . 


(21) 


Then  from  (19)  it  is  clear  that  for  P  being  P  with  V  replaced 
by  V 


V~x(y  -  X§)  =  Py, 


(22) 


so  that  the  ML  estimators  of  V  and  (3  also  can  be  re-written 
by  solving  the  following  ML  equations  simultaneously 


X'  iK1  X  (3  =  X'  V^1  y 


(23) 


and 


L  tr { V^ZjZi)  }" 


r t7/  Dt7 


(24) 


where  {c  tr(V1ZiZ/i)  }^=o  is  (N+l)xl  column  vector  with  elements 

which  are  the  trace  of  the  matrix  ( V~1ZiZli)  ,  for  i  =  0,l,...N, 
respectively. 

Also,  the  estimator  of  V  can  be  obtained  by  restricted 
maximum  likelihood  ( REML ) .  A  basic  idea  of  REML  estimation  is 
that  of  estimating  variance  components  based  on  residuals 
calculated  after  fitting  by  ordinary  least  squares  just  the 
fixed  effects  part  of  the  model.  The  REML  estimator  can  also 
by  viewed  as  maximizing  a  marginal  likelihood.  The  REML 
equations  can  therefore  be  derived  from  the  ML  equations  of 
(24)  .  By  making  suitable  replacements,  the  REML  estimator V 
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can  be  obtained  by  solving  the  following  equation: 


{  tr  [  {K'VK)  K'ZZ'K  ]  }f.0  = 


{  y'K{K'VK)  ^K'ZZ'KiK'VK)  _1 K'y  }J=0 


(25) 


with  P  =  K  (K'VK) ''  K'  where  K'X=0. 

When  the  resultant  V  is  found  by  solving  the  ML  equations 
(  (19)  ,  (20)  )  or  REML  equation  (25)  ,  the  Henderson's  mixed  model 
equation  can  be  re-written  as  follows: 


x'r'^x 

x’r'^z 

p 

X'R~xy 

z'r^x 

Z'R^Z+G'1  . 

V  . 

Z,R~1y_ 

(26) 


Finally,  p  =  [?0  ,  fj7  and  0  =  [S01,S11;  .  .  .  60W,  $1W]7  substitutes 
P  =  [Yo'Yi]7  and  v  =  [S01, 6i:l/  .  .  .  60W,  51W]/  in  (14)  as  follows: 


y  =  X  P  +  Z  v  . 


(27) 


In  fume  time  analysis,  the  test  temperature  condition  is 
constant  and  Wj  in  this  analysis  is  1.  Therefore, 
Xy=  (Wj,  Wjty)  =  ( 1 ,  ty)  and  Zjj=  ( 1 ,  ty)  .  Then  (27)  becomes 
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f  =  X§  +  Z0  =  Z0  +  Ztf 


=  Z0  +  Z((§Z/f'1(y-Zp) ) 

=  z$+zdz'v~1y-zdz,v~1z§  <28> 

=  ( J-ZGZ'V'-1)  Z0  +  {ZGZ,V~1)y 
=  (  I  -  H  )  E(  y  )  +  H  y 

where  H=ZGZ,V'1  is  the  weighted  of  mean  of  E[y]  and  y.  The 
resulting  y  is  called  a  shrinkage  estimator. 

The  predicted  value  of  fume  time  (y^)  of  lot  i  at  ty  can 

be  estimated  as 


(1,  ti:/) 


+  (1  ,ty) 


(29) 


A  90%  prediction  interval  for  a  fume  time  at  given  propellant 
age  ty  is 

y±j  ±  t  (  0.95  ,  N-2  )  SEPRED  (  |  by)  (30) 
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V'(p)  =  (  x'  V^1  X  )  -1; 


(32) 


y-(^)  =  G  Z1  V'1  G  z'  . 


In  order  to  estimate  deteriorating  patterns  of  a  group  of 
similar  lots,  one  would  be  interested  in  the  shelf-life  of  the 
group  of  similar  lots  (tsl)  that  takes  into  account  the  lot-to- 
lot  variation.  It  can  be  obtained  as  follows  by  only 
considering  the  group  effects: 


30  =  Xij  p  -  t  (  0.95  ,  N-2  )  SEPRED(gioup)  (  y  \  tsl  ) 


(33) 


=  (I/O 


-t  (0.95,  N~2 )  SEPRED 


(group) 


(y  I  taJ) 


where  SEPRED (group)  (  9  \  tsl)  is  the  standard  error  of  predicted 
y  of  a  group  given  by 


SEPRED ( gr oup ,  (  y\t±j) 


[  (l,^)7  +  t2  ]"• 


(34) 


The  group  shelf -life  can  be  compared  to  the  individual 
shelf-life  (tss,)  based  on  the  shrinkage  estimation  (28)  .  It 
is  defined  as  the  time  period  at  which  the  expected  fume  time 
curve  intersects  the  30  days  of  fume  time: 
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30  J=  Xp  +  Z<1 


(35) 


=  (  I  -  zdz'V'1'  )  zp  +  zQz'v~xy 


where  I  is  a  rijXl  vector  of  l's  and  Z  is  a  njxl  vector  of 
(l,tssl)  . 

C.  DATA  ANALYSIS 

Standard  linear  regression  model  II  based  on  the  reduced 
data  is  further  investigated  for  the  random  effects  regression 
model  analysis.  As  explained  earlier,  under  the  AAT,  the  fume 
time  ( yy )  of  each  lot  is  repeatedly  measured  five  to  nine 
times  in  an  accelerated  condition  (65.5  °C)  during  the  last  20 
years.  The  age  (ty)  for  each  of  the  lot  is  computed  by 
subtracting  the  manufacturing  date  (see  Appendix  A)  from  the 
date  in  which  when  a  45g  propellant  batch  was  sent  to  the 
heating  chamber. 

First,  in  order  to  estimate  the  safe  shelf-life  of  a 
group  of  similar  gun  propellants  based  on  the  random  effects 
regression  model,  the  within-lot  and  between-lot  models  were 
combined  to  be  a  mixed  effects  model  (14) .  Since  the  G  and  R 
(hence  V)  are  unknown  (15) ,  the  estimators  of  G  and  R  can  be 
obtained  by  a  restricted  maximum  likelihood  (REML)  (25) : 


0.07079100  -0.00450653 

-0.00450653  0.00030661. 

and  R  -  Var{e)  =  0.05926819.  Unknown  parameters  /3  and  v  can 
in  turn  be  estimated  by  Henderson's  mixed  equation  (26)  :  (1) 

P  =  =  [3.81021477,-0.14105394]’  where  se[f0]  = 

[0.07176290]  and  se [?a ]= [0 . 00534921]  ;  (2)  v  [§oi»  S117  .  .  .  §0N,  S1W] 
and  the  corresponding  standard  errors  are  shown  in  Table  9. 


G  =  Var  ( v )  = 


Qc 


'01 


°01 

0i 


30 


PARAMETER 

SUBJECT 

ESTIMATE 

STD  ERROR 

INTERCEPT 

LOT  ID 

11142 

-0 . 06196541 

0 . 19329612 

AGE 

LOT  ID 

11142 

0 . 00479259 

0 . 01346718 

INTERCEPT 

LOT  ID 

11143 

-0.11611040 

0.18778713 

AGE 

LOT  ID 

11143 

0 . 01007674 

0 . 01290832 

INTERCEPT 

LOT  ID 

11144 

-0.16769671 

0.19020538 

AGE 

LOT  ID 

11144 

0 . 01415533 

0 . 01316629 

INTERCEPT 

LOT  ID 

11145 

-0 . 09516341 

0.19031856 

AGE 

LOT  ID 

11145 

0 . 00494655 

0.01335297 

INTERCEPT 

LOT  ID 

11146 

-0 . 08786149 

0 . 18872403 

AGE 

LOT  ID 

11146 

0 . 00599071 

0 . 01323274 

INTERCEPT 

LOT  ID 

11147 

-0.21248452 

0.18231704 

AGE 

LOT  ID 

11147 

0 . 01167817 

0 . 01275956 

INTERCEPT 

LOT  ID 

11148 

-0.11397421 

0.18758544 

AGE 

LOT  ID 

11148 

0 . 00853627 

0 . 01298120 

INTERCEPT 

LOT  ID 

11149 

-0 . 07768438 

0.18631350 

AGE 

LOT  ID 

11149 

-0 . 00123725 

0 . 01287421 

INTERCEPT 

LOT  ID 

11150 

-0 . 15424392 

0.18762834 

AGE 

LOT  ID 

11150 

0 . 00950736 

0 . 01304139 

INTERCEPT 

LOT  ID 

11151 

-0.12749461 

0.18679008 

AGE 

LOT  ID 

11151 

0 . 00861269 

0 . 01295106 

INTERCEPT 

LOT  ID 

11156 

-0.15799910 

0.18658553 

AGE 

LOT  ID 

11156 

0 . 00818934 

0 . 01299146 

INTERCEPT 

LOT  ID 

11157 

-0 . 17588616 

0.18545442 

AGE 

LOT  ID 

11157 

0 . 01418657 

0.01281651 

INTERCEPT 

LOT  ID 

11158 

-0.13427054 

0 . 18506946 

AGE 

LOT  ID 

11158 

0 . 00509910 

0 . 01293535 

INTERCEPT 

LOT  ID 

11159 

-0.23813672 

0.18472938 

AGE 

LOT  ID 

11159 

0 . 01251881 

0 . 01296406 

INTERCEPT 

LOT  ID 

11160 

-0.22564272 

0.18029667 

AGE 

LOTID 

11160 

0 . 01415942 

0 . 01265016 

INTERCEPT 

LOT  ID 

11181 

-0.16269825 

0.18241499 

AGE 

LOTID 

11181 

0 . 00759863 

0 . 01279735 

INTERCEPT 

LOTID 

11182 

-0 . 14882165 

0 . 18459049 

AGE 

LOTID 

11182 

0 . 01126227 

0 . 01279876 

INTERCEPT 

LOTID 

11183 

0.17988730 

0.18660005 

AGE 

LOTID 

11183 

-0 . 01074349 

0 . 01305771 

INTERCEPT 

LOTID 

11184 

0.21035950 

0.18618665 

AGE 

LOTID 

11184 

-0 . 01320753 

0 . 01303212 

INTERCEPT 

LOTID 

11185 

0.44493213 

0.19631611 

AGE 

LOTID 

11185 

-0 . 02624763 

0 . 01352684 

INTERCEPT 

LOTID 

11188 

0.24783741 

0 . 18929810 

AGE 

LOTID 

11188 

-0 . 01470002 

0 . 01319538 

INTERCEPT 

LOTID 

11189 

0.22284679 

0.18633335 

AGE 

LOTID 

11189 

-0 . 01325530 

0 . 01303528 

INTERCEPT 

LOTID 

11190 

0 . 18150475 

0 . 18692804 

AGE 

LOTID 

11190 

-0 . 01253802 

0 . 01310098 

INTERCEPT 

LOTID 

11191 

0.22307183 

0 . 18628924 

AGE 

LOTID 

11191 

-0 . 01355682 

0 . 01303573 

INTERCEPT 

LOTID 

11192 

0.20898367 

0 . 18614745 
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AGE 

INTERCEPT 

AGE 

INTERCEPT 

AGE 

INTERCEPT 
AGE _ 


LOTID  11192 
LOTID  11193 
LOTID  11193 
LOTID  11194 
LOTID  11194 
LOTID  11195 
LOTID  11195 


-0 . 01170315 
0.19615700 
-0 . 01181337 
0.17762001 
-0.01152574 
0.16493381 
-0 . 01078221 


0.01302592 
0.18504953 
0.01297891 
0.18452648 
0.01310043 
0.18380045 
0 . 01291217 


Table  9 .  $  And  Std  Error  of  $ 


Subsequently,  the  safe  shelf-life  of  the  group  of  similar  lots 
(tsl)  that  takes  into  account  the  lot-to-lot  variation  can  be 
obtained  as  follows: 


30= (1, 


./  3.81021477 
sJ  \  -0.14105394 


-C(0.9  5,N-2)SEPRED{group] 


(y 


Figure  C.22  shows  the  resulting  shelf -life  is  8405.12  days 
(23.02  years) . 

This  result  can  be  compared  to  the  one  obtain  by  the 
current  method  (5) ,  which  is  8469  days  (23.20  years) .  Note 
that  the  current  method  without  taking  into  account  the 
potential  lot-to-lot  variation  provides  a  liberal  shelf-life 
that  is  64  days  longer  than  the  random  effects  regression 
model . 

Finally,  individual  shelf-life  of  lot  i  (tssl)  based  on  the 
shrinkage  estimation  (28)  is  computed  by  checking  the  time 
period  at  which  the  expected  fume  time  curve  intersects  the  30 
days  of  fume  time  (see  Table  10)  .  The  resulting  individual 
shelf -life  (tss,)  varies  from  27.998  to  24.942  years.  They  all 
exceed  the  group  shelf -life  based  on  either  the  fixed  effects 
(23.20  years)  or  the  random  effects  (23.02  years)  regression 
model . 

The  group  shelf-lives  based  on  the  prediction  limit  for 
the  fume  time  can  be  compared  to  those  based  on  the  confidence 
limit  for  the  expected  value.  They  turn  out  to  be  25.66  and 
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OBS  LOTID  shelf-life  shelf-life  shelf-life 
(shrinkage  (confidence  (prediction 
estimation)  limit)  limit) 


1 

11142 

26 . 914 

24.520 

23.427 

2 

11143 

27.576 

25.131 

23 . 971 

3 

11144 

28 . 056 

25.411 

24.315 

4 

11145 

26 . 691 

24.365 

23.219 

5 

11146 

26 . 951 

24.726 

23.424 

6 

11147 

27.173 

24 . 729 

23.561 

7 

11148 

27 . 272 

24 . 863 

23 . 698 

8 

11149 

25.654 

23 . 630 

22.459 

9 

11150 

27.167 

24.726 

23.561 

10 

11151 

27.185 

24.794 

23.629 

11 

11156 

27 . 074 

24.589 

23.356 

12 

11157 

27 . 998 

25.411 

24.246 

13 

11158 

26.433 

24 . 178 

23 . 013 

14 

11159 

27 . 151 

24.522 

23.492 

15 

11160 

27 . 600 

25 . 069 

23 . 904 

16 

11181 

26 . 715 

24.383 

23.218 

17 

11182 

27.576 

25 . 068 

23 . 972 

18 

11183 

25.563 

23.767 

22 . 663 

19 

11184 

25.534 

23 . 629 

22.552 

20 

11185 

24 . 942 

23.212 

22 . 191 

21 

11188 

25.526 

23 . 631 

22.573 

22 

11189 

25 . 603 

23 . 698 

22.621 

23 

11190 

25.275 

23.493 

22.459 

24 

11191 

25.555 

23 . 632 

22.578 

25 

11192 

25 . 772 

23 . 835 

22 . 735 

26 

11193 

25 . 672 

23 . 771 

22 . 658 

27 

11194 

25 . 590 

23.561 

22.550 

28 

11195 

25.639 

23.698 

22 . 614 

Table  10.  Individual  Shelf -Life  of  Propellant  Lot 

24.314  years,  which  are  longer  than  their  counterparts  based 
on  the  prediction  limit. 

The  individual  shelf-lives  can  also  be  calculated  by 
checking  the  time  period  when  95%  one-sided  lower  prediction 
or  confidence  limits  intersect  30  days  of  fume  time.  The 
results  are  shown  in  Table  10.  About  42%  of  individual  shelf- 
lives  based  on  the  prediction  limit  turned  out  to  be  shorter 
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than  the  group  shelf-life  (23.02  years)  .  The  rest  are  as  same 
as  the  shrinkage  estimation  that  all  exceeds  the  group  shelf- 
life  based  on  either  the  fixed  and  random  effects  regression 
model . 

All  calculations  are  done  by  using  PROC  MIXED  and  PROC 
REG  of  a  statistical  package  SAS  [Ref. 12]. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Most  gun  propellant  is  stored  at  depots  for  a  long  time 
before  it  is  used.  In  storage,  the  quality  of  the  gun 
propellant  may  deteriorate  and  become  unstable.  In  an  attempt 
to  avoid  disaster  due  to  use  of  unstable  gun  propellant, 
accurate  prediction  of  the  safe  shelf-life  of  gun  propellant 
is  necessary.  The  shelf-life  estimation  methods  used 
currently  for  a  group  of  similar  gun  propellant  lots  are  based 
on  a  fixed  effects  regression  model.  This  does  not  take  into 
consideration  the  fact  that  samples  from  the  same  lot  are  more 
similar  than  samples  between  lots.  To  capitalize  on  this  lot¬ 
to-lot  variation  when  estimating  the  shelf-life,  first,  a 
random  effects  regression  model  (RERM)  is  developed. 
Secondly,  a  combined  mixed  effects  model  is  estimated.  The 
estimated  model  is  then  used  to  predict  not  only  the  shelf- 
life  of  a  group  of  similar  lots  but  also  that  of  each 
individual  lot  of  5"/54  NACO  gun  propellant  stockpile. 
Finally,  results  of  the  data  analysis  are  discussed  and 
recommendations  are  made. 

The  claimed  shelf-life  of  the  NACO  propellant  lots 
provided  by  NSWC  is  35  years.  It  is  interesting  to  note  that 
when  the  confidence  limit  method  is  used  instead  of  the 
prediction  limit  method,  the  standard  linear  regression  model 

I  provides  35  years  as  the  estimated  shelf -life.  However,  the 
assumptions  for  standard  linear  regression  model  I  were  not 
validated.  Therefore,  the  use  of  the  claimed  shelf-life  (35 
years)  is  not  supported  and  needs  to  be  amended. 

The  group  shelf- life  estimated  based  on  a  random  effects 
regression  model  turns  out  to  be  23.02  years.  A  fixed  effects 
regression  model  based  on  the  standard  linear  regression  model 

II  provides  a  slightly  longer  shelf-life  (23.20  years)  than 
23.02  years.  The  individual  shelf -lives  based  on  the 


35 


shrinkage  estimation  vary  from  27.998  to  24.942  years  which 
all  exceed  the  group  shelf -lives  obtained  from  both  the  fixed 
and  the  random  effects  regression  models.  It  indicates  that 
group  shelf-life  estimated  based  on  the  fixed  or  the  random 
effects  regression  models  can  be  relatively  conservative 
compared  to  the  individual  shelf -lives.  When  the  group  shelf - 
life  (23.02  years)  is  applied  to  the  maintenance  of  unstable 
stockpiles,  there  is  a  chance  that  all  individual  lots  would 
be  discarded  even  though  they  are  still  safe.  In  view  of 
such  opportunity  loss,  use  of  individual  shelf -life  is 
recommended  over  the  group  shelf-life  when  the  management  of 
individual  lot  is  expected  to  cost  less  than  that  of  the  bulk 
management . 

B .  RECOMMENDATIONS 

In  this  thesis,  a  particular  propellant  stockpile  data 
set  is  analyzed.  Although  the  approaches  proposed  can  be 
applied  to  any  other  types  of  propellant,  the  model  selection 
may  be  dictated  by  the  nature  of  data  set. 

There  is  a  nonlinear  decreasing  trend  in  fume  time  in 
some  individual  lots  as  the  propellant  ages.  This  raises 
questions  about  the  approaches  that  use  traditional  linear 
least  squares  analysis  of  fume  time  data  to  predict  safe 
shelf -life.  This  thesis  recommends  that  the  nonlinear 
regression  analysis  approach  be  conducted  for  some  individual 
propellant  lots. 
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APPENDIX  A  MANUFACTURING  DATE 


OBS 

LOTID 

MANUFACTURING 

OBS 

LOTID 

MANUFACTURING 

DATE 

DATE 

1 

11142 

10/16/1972 

2 

11143 

11/06/1972 

3 

11144 

11/07/1972 

4 

11145 

11/28/1972 

5 

11146 

12/07/1972 

6 

11147 

4/01/1973 

7 

11148 

12/27/1972 

8 

11149 

12/27/1972 

9 

11150 

12/27/1972 

10 

11151 

1/24/1973 

11 

11156 

1/24/1973 

12 

11157 

2/15/1973 

13 

11158 

4/01/1973 

14 

11159 

4/01/1973 

15 

11160 

10/18/1973 

16 

11181 

10/19/1973 

17 

11182 

4/01/1973 

18 

11183 

10/24/1973 

19 

11184 

10/25/1973 

20 

11185 

7/26/1973 

21 

11188 

7/29/1973 

22 

11189 

10/31/1973 

23 

11190 

10/31/1973 

24 

11191 

11/01/1973 

25 

11192 

11/06/1973 

26 

11193 

1/14/1974 

27 

11194 

1/14/1974 

28 

11195 

1/14/1974 
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APPENDIX  B  COMPUTER  PROGRAM  LIST 


The  following  computer  program  is  a  SAS  code. 


OPTIONS  LS 

=  80/ 

/* 

MASTER  SAMPLES 

*/ 

DATA  ALL1 ; 

INPUT 

OB 

ID 

INDEX 

EAGE 

EDTF 

EYAGE 

EYDTF ; 

ID=  N  ; 

CARDS ; 

1 

1 

11142 

204 

756 

0.5589 

2.07123 

2 

1 

11142 

960 

875 

2.6301 

2.39726 

3 

1 

11142 

1834 

1117 

5 . 0247 

3 . 06027 

4 

1 

11142 

2951 

1002 

8 . 0849 

2.74521 

5 

1 

11142 

3952 

741 

10.8274 

2.03014 

6 

1 

11142 

4692 

688 

12.8548 

1.88493 

7 

1 

11142 

5380 

710 

14 . 7397 

1.94521 

8 

1 

11142 

6238 

578 

17.0904 

1.58356 

9 

2 

11143 

184 

724 

0 . 5041 

1.98356 

10 

2 

11143 

908 

837 

2 .4877 

2.29315 

11 

2 

11143 

1744 

1083 

4.7781 

2 . 96712 

12 

2 

11143 

2827 

1058 

7 . 7452 

2 .89863 

13 

2 

11143 

3884 

761 

10 . 6411 

2 . 08493 

14 

2 

11143 

4644 

699 

12 . 7233 

1.91507 

15 

2 

11143 

5343 

727 

14.6384 

1.99178 

16 

2 

11143 

6218 

589 

17.0356 

1.61370 

17 

2 

11143 

6807 

552 

18.6493 

1.51233 

18 

3 

11144 

276 

663 

0.7562 

1.81644 

19 

3 

11144 

939 

833 

2.5726 

2.28219 

20 

3 

11144 

1771 

1077 

4 . 8521 

2 . 95068 

21 

3 

11144 

2848 

1043 

7.8027 

2.85753 

22 

3 

11144 

3890 

742 

10.6575 

2.03288 

23 

3 

11144 

4631 

683 

12.6877 

1.87123 

24 

3 

11144 

5314 

685 

14.5589 

1.87671 

25 

3 

11144 

5998 

552 

16.4329 

1 . 51233 

26 

3 

11144 

6363 

806 

17.4329 

2.20822 

27 

4 

11145 

913 

808 

2.5014 

2.21370 

28 

4 

11145 

1720 

1077 

4.7123 

2.95068 

29 

4 

11145 

2797 

1046 

7 . 6630 

2 .86575 

30 

4 

11145 

3842 

764 

10.5260 

2 .09315 

31 

4 

11145 

4605 

683 

12 . 6164 

1.87123 

32 

4 

11145 

5288 

613 

14.4877 

1.67945 

33 

4 

11145 

5900 

559 

16 . 1644 

1.53151 

34 

5 

11146 

152 

689 

0.4164 

1.88767 

35 

5 

11146 

841 

834 

2.3041 

2.28493 

36 

5 

11146 

1674 

1045 

4.5863 

2.86301 

37 

5 

11146 

2719 

1106 

7.4493 

3 . 03014 

38 

5 

11146 

3824 

111 

10.4767 

2.12877 

39 

5 

11146 

4600 

694 

12 .6027 

1.90137 

40 

5 

11146 

5294 

717 

14.5041 

1.96438 

41 

5 

11146 

6010 

556 

16.4658 

1.52329 

42 

6 

11147 

0 

661 

0 . 0000 

1.81096 

43 

6 

11147 

661 

855 

1 . 8110 

2.34247 

44 

6 

11147 

1515 

1017 

4.1507 

2 . 78630 

45 

6 

11147 

2532 

1057 

6 . 9370 

2 .89589 

46 

6 

11147 

3588 

786 

9.8301 

2.15342 

47 

6 

11147 

4373 

700 

11.9808 

1.91781 

39 


48 

6 

11147 

5073 

655 

13 . 8986 

1.79452 

49 

6 

11147 

5727 

560 

15.6904 

1.53425 

50 

6 

11147 

6287 

459 

17.2247 

1.25753 

51 

7 

11148 

132 

728 

0.3616 

1.99452 

52 

7 

11148 

857 

834 

2.3479 

2.28493 

53 

7 

11148 

1690 

1077 

4.6301 

2 . 95068 

54 

7 

11148 

2767 

1065 

7.5808 

2 . 91781 

55 

7 

11148 

3831 

761 

10.4959 

2.08493 

56 

7 

11148 

4591 

708 

12 . 5781 

1.93973 

57 

7 

11148 

5299 

696 

14.5178 

1.90685 

58 

7 

11148 

5994 

528 

16.4219 

1.44658 

59 

7 

11148 

6522 

596 

17 . 8685 

1.63288 

60 

8 

11149 

888 

806 

2.4329 

2.20822 

61 

8 

11149 

1693 

1046 

4.6384 

2.86575 

62 

8 

11149 

2739 

1086 

7.5041 

2.97534 

63 

8 

11149 

3824 

737 

10.4767 

2 . 01918 

64 

8 

11149 

4560 

699 

12.4932 

1.91507 

65 

8 

11149 

5259 

637 

14.4082 

1.74521 

66 

8 

11149 

5895 

19 

16 . 1507 

0.05205 

67 

8 

11149 

6813 

477 

18.6658 

1.30685 

68 

9 

11150 

132 

728 

0 . 3616 

1.99452 

69 

9 

11150 

860 

824 

2.3562 

2.25753 

70 

9 

11150 

1683 

1063 

4.6110 

2 . 91233 

71 

9 

11150 

2746 

1074 

7 . 5233 

2.94247 

72 

9 

11150 

3819 

728 

10.4630 

1.99452 

73 

9 

11150 

4546 

682 

12.4548 

1.86849 

74 

9 

11150 

5228 

512 

14.3233 

1.40274 

75 

9 

11150 

5740 

702 

15 . 7260 

1 . 92329 

76 

9 

11150 

6441 

538 

17.6466 

1.47397 

77 

10 

11151 

844 

823 

2.3123 

2.25479 

78 

10 

11151 

1666 

1056 

4.5644 

2.89315 

79 

10 

11151 

2722 

1083 

7.4575 

2 . 96712 

80 

10 

11151 

3804 

747 

10.4219 

2 . 04658 

81 

10 

11151 

4550 

698 

12.4658 

1.91233 

82 

10 

11151 

5248 

711 

14.3781 

1.94795 

83 

10 

11151 

5958 

523 

16.3233 

1.43288 

84 

10 

11151 

6481 

557 

17.7562 

1.52603 

85 

11 

11156 

106 

709 

0.2904 

1.94247 

86 

11 

11156 

815 

842 

2.2329 

2.30685 

87 

11 

11156 

1656 

1038 

4.5370 

2 . 84384 

88 

11 

11156 

2694 

1062 

7.3808 

2 . 90959 

89 

11 

11156 

3755 

765 

10.2877 

2.09589 

90 

11 

11156 

4519 

665 

12.3808 

1.82192 

91 

11 

11156 

5184 

589 

14.2027 

1.61370 

92 

11 

11156 

5773 

592 

15 . 8164 

1.62192 

93 

11 

11156 

6364 

463 

17.4356 

1.26849 

94 

12 

11157 

822 

822 

2.2521 

2.25205 

95 

12 

11157 

1643 

1084 

4 . 5014 

2 . 96986 

96 

12 

11157 

2727 

1056 

7.4712 

2.89315 

97 

12 

11157 

3782 

746 

10.3616 

2 . 04384 

98 

12 

11157 

4527 

710 

12.4027 

1 . 94521 

99 

12 

11157 

5237 

723 

14.3479 

1.98082 

100 

12 

11157 

6118 

563 

16 . 7616 

1 . 54247 

101 

12 

11157 

6681 

662 

18 .3041 

1 . 81370 

102 

13 

11158 

0 

738 

0.0000 

2 . 02192 

103 

13 

11158 

738 

829 

2 . 0219 

2.27123 

104 

13 

11158 

1566 

1092 

4.2904 

2.99178 

105 

13 

11158 

2658 

1041 

7.2822 

2.85205 

106 

13 

11158 

3698 

746 

10 . 1315 

2 . 04384 

107 

13 

11158 

4443 

681 

12 . 1726 

1.86575 

108 

13 

11158 

5124 

620 

14 . 0384 

1.69863 

40 


109 

13 

11158 

5743 

448 

15.7342 

1.22740 

110 

13 

11158 

6191 

482 

16.9616 

1.32055 

111 

14 

11159 

0 

730 

0.0000 

2 . 00000 

112 

14 

11159 

730 

842 

2 . 0000 

2.30685 

113 

14 

11159 

1571 

1116 

4.3041 

3 . 05753 

114 

14 

11159 

2687 

1042 

7.3616 

2 . 85479 

115 

14 

11159 

3728 

550 

10.2137 

1.50685 

116 

14 

11159 

4277 

624 

11.7178 

1.70959 

117 

14 

11159 

4901 

661 

13.4274 

1.81096 

118 

14 

11159 

5562 

666 

15.2384 

1.82466 

119 

14 

11159 

6227 

489 

17.0603 

1.33973 

120 

15 

11160 

582 

834 

1.5945 

2.28493 

121 

15 

11160 

1415 

1085 

3 . 8767 

2 . 97260 

122 

15 

11160 

2500 

1038 

6 . 8493 

2 . 84384 

123 

15 

11160 

3537 

745 

9.6904 

2 . 04110 

124 

15 

11160 

4281 

695 

11.7288 

1 . 90411 

125 

15 

11160 

4976 

720 

13.6329 

1.97260 

126 

15 

11160 

5695 

573 

15.6027 

1.56986 

127 

15 

11160 

6268 

581 

17.1726 

1.59178 

128 

16 

11181 

632 

826 

1.7315 

2.26301 

129 

16 

11181 

1457 

1177 

3 . 9918 

3.22466 

130 

16 

11181 

2633 

953 

7.2137 

2.61096 

131 

16 

11181 

3586 

721 

9 . 8247 

1 . 97534 

132 

16 

11181 

4306 

680 

11 .7973 

1 . 86301 

133 

16 

11181 

4986 

686 

13.6603 

1 . 87945 

134 

16 

11181 

5671 

524 

15.5370 

1.43562 

135 

16 

11181 

6195 

488 

16 . 9726 

1.33699 

136 

17 

11182 

0 

741 

0 . 0000 

2 . 03014 

137 

17 

11182 

741 

841 

2 . 0301 

2.30411 

138 

17 

11182 

1581 

1131 

4.3315 

3 . 09863 

139 

17 

11182 

2712 

1023 

7.4301 

2.80274 

140 

17 

11182 

3734 

747 

10.2301 

2.04658 

141 

17 

11182 

4480 

691 

12.2740 

1.89315 

142 

17 

11182 

5171 

716 

14.1671 

1.96164 

143 

17 

11182 

6008 

595 

16 .4603 

1.63014 

144 

17 

11182 

6603 

600 

18 . 0904 

1 . 64384 

145 

18 

11183 

622 

840 

1 . 7041 

2.30137 

146 

18 

11183 

1461 

1300 

4 . 0027 

3 . 56164 

147 

18 

11183 

2760 

1105 

7.5616 

3 . 02740 

148 

18 

11183 

3865 

779 

10 . 5890 

2.13425 

149 

18 

11183 

4643 

734 

12 . 7205 

2 . 01096 

150 

18 

11183 

5377 

713 

14 . 7315 

1.95342 

151 

18 

11183 

6089 

539 

16 . 6822 

1.47671 

152 

19 

11184 

586 

846 

1.6055 

2.31781 

153 

19 

11184 

1431 

1322 

3 . 9205 

3.62192 

154 

19 

11184 

2752 

1112 

7 . 5397 

3 . 04658 

155 

19 

11184 

3864 

781 

10.5863 

2.13973 

156 

19 

11184 

4644 

736 

12 . 7233 

2 . 01644 

157 

19 

11184 

5380 

692 

14 . 7397 

1.89589 

158 

19 

11184 

6071 

515 

16.6329 

1.41096 

159 

20 

11185 

0 

698 

0 . 0000 

1.91233 

160 

20 

11185 

698 

1049 

1.9123 

2 .87397 

161 

20 

11185 

1746 

1492 

4 . 7836 

4 . 08767 

162 

20 

11185 

3237 

998 

8 . 8685 

2 . 73425 

163 

20 

11185 

4234 

784 

11 .6000 

2 . 14795 

164 

20 

11185 

5018 

766 

13 . 7479 

2 . 09863 

165 

20 

11185 

5783 

569 

15.8438 

1.55890 

166 

20 

11185 

6352 

579 

17 .4027 

1 . 58630 

167 

21 

11188 

0 

684 

0.0000 

1.87397 

168 

21 

11188 

684 

854 

1 . 8740 

2.33973 

169 

21 

11188 

1537 

1358 

4.2110 

3 . 72055 

41 


170 

21 

11188 

2894 

1072 

7.9288 

2.93699 

171 

21 

11188 

3966 

766 

10.8658 

2.09863 

172 

21 

11188 

4731 

743 

12.9616 

2.03562 

173 

21 

11188 

5474 

689 

14 . 9973 

1.88767 

174 

21 

11188 

6162 

553 

16.8822 

1.51507 

175 

22 

11189 

597 

831 

1.6356 

2.27671 

176 

22 

11189 

1427 

1351 

3 . 9096 

3 . 70137 

177 

22 

11189 

2777 

1085 

7.6082 

2.97260 

178 

22 

11189 

3862 

795 

10.5808 

2 . 17808 

179 

22 

11189 

4656 

740 

12.7562 

2.02740 

180 

22 

11189 

5396 

695 

14.7836 

1.90411 

181 

22 

11189 

6090 

564 

16.6849 

1.54521 

182 

23 

11190 

596 

866 

1.6329 

2.37260 

183 

23 

11190 

1461 

1316 

4 . 0027 

3 . 60548 

184 

23 

11190 

2776 

1075 

7.6055 

2.94521 

185 

23 

11190 

3851 

798 

10.5507 

2.18630 

186 

23 

11190 

4648 

719 

12 . 7342 

1.96986 

187 

23 

11190 

5367 

586 

14 . 7041 

1.60548 

188 

23 

11190 

5952 

540 

16.3068 

1.47945 

189 

24 

11191 

601 

825 

1.6466 

2.26027 

190 

24 

11191 

1425 

1349 

3 . 9041 

3.69589 

191 

24 

11191 

2773 

1096 

7.5973 

3 . 00274 

192 

24 

11191 

3869 

770 

10.6000 

2.10959 

193 

24 

11191 

4638 

757 

12 . 7068 

2 . 07397 

194 

24 

11191 

5395 

683 

14 . 7808 

1.87123 

195 

24 

11191 

6077 

553 

16 . 6493 

1.51507 

196 

25 

11192 

598 

824 

1.6384 

2.25753 

197 

25 

11192 

1421 

1357 

3 . 8932 

3 . 71781 

198 

25 

11192 

2777 

1069 

7.6082 

2 . 92877 

199 

25 

11192 

3846 

799 

10.5370 

2 . 18904 

200 

25 

11192 

4644 

746 

12.7233 

2 . 04384 

201 

25 

11192 

5390 

700 

14.7671 

1.91781 

202 

25 

11192 

6089 

615 

16 .6822 

1.68493 

203 

26 

11193 

528 

851 

1.4466 

2 . 33151 

204 

26 

11193 

1378 

1361 

3.7753 

3 . 72877 

205 

26 

11193 

2738 

1076 

7 . 5014 

2 . 94795 

206 

26 

11193 

3814 

764 

10.4493 

2 . 09315 

207 

26 

11193 

4577 

749 

12.5397 

2 . 05205 

208 

26 

11193 

5326 

699 

14.5918 

1.91507 

209 

26 

11193 

6024 

585 

16 . 5041 

1.60274 

210 

27 

11194 

528 

825 

1.4466 

2.26027 

211 

27 

11194 

1352 

1329 

3 . 7041 

3.64110 

212 

27 

11194 

2680 

1075 

7.3425 

2 . 94521 

213 

27 

11194 

3755 

828 

10.2877 

2.26849 

214 

27 

11194 

4582 

742 

12.5534 

2 . 03288 

215 

28 

11195 

528 

823 

1.4466 

2.25479 

216 

28 

11195 

1350 

1331 

3.6986 

3.64658 

217 

28 

11195 

2680 

1059 

7.3425 

2 . 90137 

218 

28 

11195 

3739 

808 

10.2438 

2.21370 

219 

28 

11195 

4546 

753 

12.4548 

2 . 06301 

220 

28 

11195 

5299 

694 

14.5178 

1 . 90137 

221 

28 

11195 

5992 

507 

16.4164 

1.38904 

/* 

FLEET  RETURN 

DATA 

*/ 

DATA  ALL; 

INPUT 

OBS 

ID 

INDEX 

MDATE 

DTF 

AGE 

CARD; 

1 

1 

11142 

72.8795 

716 

1671.0 

2 

1 

11142 

72 . 8795 

662 

4372 . 0 

3 

1 

11142 

72 . 8795 

812 

1671 . 0 

42 


65 

7 

11148 

73.0767 

707 

2801.0 

66 

7 

11148 

73.0767 

809 

1599.0 

67 

7 

11148 

73.0767 

681 

2433.0 

68 

7 

11148 

73.0767 

645 

4300.0 

69 

7 

11148 

73.0767 

703 

2824.0 

70 

8 

11149 

73 . 0767 

791 

1599.0 

71 

8 

11149 

73 . 0767 

628 

4300.0 

72 

8 

11149 

73 . 0767 

638 

4300.0 

73 

8 

11149 

73.0767 

611 

5004.0 

74 

8 

11149 

73.0767 

606 

5045.0 

75 

9 

11150 

73 . 0767 

667 

2801.0 

76 

9 

11150 

73 . 0767 

673 

2801.0 

77 

9 

11150 

73 . 0767 

686 

2801.0 

78 

9 

11150 

73 . 0767 

606 

5045.0 

79 

9 

11150 

73.0767 

594 

2234.0 

80 

9 

11150 

73.0767 

628 

4300.0 

81 

9 

11150 

73.0767 

634 

2103 . 0 

82 

9 

11150 

73.0767 

809 

1599.0 

83 

9 

11150 

73 . 0767 

635 

4300.0 

84 

10 

11151 

73.1493 

612 

4977.5 

85 

10 

11151 

73.1493 

611 

4290.5 

86 

10 

11151 

73.1493 

638 

4273.5 

87 

10 

11151 

73.1493 

896 

751.5 

88 

10 

11151 

73 . 1493 

760 

2650 . 5 

89 

10 

11151 

73.1493 

797 

1444.5 

90 

10 

11151 

73.1493 

766 

2671.5 

91 

10 

11151 

73.1493 

615 

5018 . 5 

92 

10 

11151 

73.1493 

860 

943 . 5 

93 

10 

11151 

73 . 1493 

692 

2523 . 5 

94 

11 

11156 

73.1493 

678 

4160.5 

95 

11 

11156 

73 . 1493 

684 

4160 . 5 

96 

11 

11156 

73.1493 

590 

5018.5 

97 

11 

11156 

73.1493 

604 

4290.5 

98 

11 

11156 

73 . 1493 

697 

1302 . 5 

99 

11 

11156 

73 . 1493 

615 

5018.5 

100 

12 

11157 

73.2082 

613 

4997.0 

101 

12 

11157 

73.2082 

610 

4956 . 0 

102 

12 

11157 

73.2082 

623 

4252 . 0 

103 

12 

11157 

73.2082 

782 

1907 . 0 

104 

12 

11157 

73.2082 

638 

4252.0 

105 

12 

11157 

73 .2082 

627 

2055.0 

106 

12 

11157 

73.2082 

667 

2378 . 0 

107 

12 

11157 

73.2082 

872 

971.0 

108 

12 

11157 

73.2082 

806 

681.0 

109 

13 

11158 

73.3370 

624 

4909.0 

110 

13 

11158 

73 . 3370 

697 

2664.0 

111 

13 

11158 

73.3370 

693 

2686 . 0 

112 

13 

11158 

73.3370 

707 

2686 . 0 

113 

13 

11158 

73 . 3370 

627 

4222 . 0 

114 

13 

11158 

73.3370 

675 

2338 . 0 

115 

13 

11158 

73 . 3370 

649 

4205.0 

116 

13 

11158 

73.3370 

728 

1635.0 

117 

13 

11158 

73 .3370 

615 

4950.0 

118 

14 

11159 

73 . 3370 

612 

4909 . 0 

119 

14 

11159 

73.3370 

616 

4950.0 

120 

14 

11159 

73 .3370 

930 

924 . 0 

121 

15 

11160 

73.8849 

744 

434.0 

122 

15 

11160 

73 . 8849 

728 

1435 . 0 

123 

15 

11160 

73 . 8849 

672 

1082.0 

124 

15 

11160 

73 . 8849 

677 

4005.0 

125 

15 

11160 

73 . 8849 

690 

4005 . 0 

44 
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187 

26 

11193 

74.1219 

702 

2399.5 

188 

26 

11193 

74 . 1219 

703 

2399.5 

189 

26 

11193 

74.1219 

710 

2399.5 

190 

26 

11193 

74.1219 

726 

2399.5 

191 

26 

11193 

74 . 1219 

918 

637.5 

192 

26 

11193 

74.1219 

711 

2399.5 

193 

26 

11193 

74.1219 

725 

2399.5 

194 

26 

11193 

74.1219 

726 

2399.5 

195 

26 

11193 

74.1219 

731 

2399.5 

196 

26 

11193 

74.1219 

697 

2044.5 

197 

26 

11193 

74.1219 

700 

2051.5 

198 

26 

11193 

74.1219 

716 

1217.5 

199 

26 

11193 

74 . 1219 

654 

3918 . 5 

200 

26 

11193 

74 . 1219 

646 

3918.5 

201 

27 

11194 

74.1219 

647 

4622 . 5 

202 

27 

11194 

74.1219 

688 

2044.5 

203 

28 

11195 

74.1219 

612 

4622.5 

204 

28 

11195 

74 . 1219 

609 

4622.5 

PROC  SORT;  BY  INDEX; 

PROC  PRINT  DATA=  ALL1 ; VAR  INDEX  EAGE  EDTF  EYAGE  EYDTF ; 

/*  COMBINE  MASTER  SAMPLE  AND  FLEET  RETURN  DATA  */ 

DATA  NEW; 

MERGE  ALL1  ALL 2  ;  BY  INDEX; 

IF  ID  =  LAG (ID)  THEN  EAGE  = 

IF  ID  =  LAG (ID)  THEN  EDTF  =  .; 

PROC  PRINT  DATA=NEW;  VAR  INDEX  MDATE  DTF  AGE  EDTF  EAGE; 

PROC  UNIVARIATE; 

PROC  PLOT; 

PLOT  DTF*AGE= ' A'  EDTF*EAGE= '  *  '  ; 

PROC  PLOT; 

PLOT  DTF*AGE= ' A'  EDTF*EAGE= ' * ' /OVERLAY; 

/*  STANDARD  LINEAR  REGRESSION  MODEL  I  */ 

PROC  REG  ;  MODEL  DTF  =  AGE  ; 

/*  STANDARD  LINEAR  REGRESSION  MODEL  II  */ 

/*  (  DELET  EAGE  BELOW  1100  DAYS  )  */ 


DATA  ALL3 (DROP=  OB  ID  EYAGE  EYDTF); SET  ALL1 ; 

DATA  ALL4 ; SET  ALL 3 ; 

IF  EAGE  LE  1100  THEN  DELETE; 

PROC  PLOT  ; 

PLOT  EDTF*EAGE= '  *  '  ; 

PROC  REG  ;  MODEL  EDTF=  EAGE; 

OUTPUT  OUT = OUT 1  P=  PRED  R=RES ; 

PROC  PRINT; 

DATA  ALL5  (DROP=  RES  )  ; MERGE  ALL4  OUT1 ;  BY  INDEX; 
PROC  PRINT; 

DATA  ALL6 (DROP=OBS) ;  SET  ALL 5 ;BY  INDEX; 

W= (EAGE-4146 . 15) * (EAGE-4146 .15) ; 

SXX=  2521553* (180-1) ; 

SEPRED  =  98.295  *  SQRT (1+ (1/180) + (W/SXX) ) ; 
L95=PRED- 1 . 645*SEPRED; 
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U95=PRED+1 . 645*SEPRED; 

EAGE30=30; 

PROC  PRINT; 

PROC  PLOT  ; 

PLOT  EDTF*EAGE= ' A '  PRED*EAGE= ' - '  L95*EAGE= ' . ' 

U95*EAGE= '  . '  EAGE3  0  *EAGE= ' - ' /OVERLAY; 

PLOT  PRED*EAGE= ' * '  L95*EAGE= ' . ' 

U95*EAGE= '  . '  EAGE 3  0  * E AGE  = ' - ' /OVERLAY; 

DATA  EXTEN1 ; 

DO  DAGE  =  8000  TO  12000  BY  100; 

PRED=1045 . 510-0 . 07*EAGE; 

W= (EAGE -4 14 6 . 15) * (EAGE-4146 .  15)  ; 

SXX=  2521553*  (180-1)  ; 

SEPRED  =  98.295  *  SQRT (1+ (1/180) + (W/SXX) ) ; 

L95=PRED- 1 . 645*SEPRED; 

U95=PRED+1 . 645*SEPRED; 

OUTPUT ; 

END; 

PROC  PRINT  DATA=EXTEN1 ; VAR  EAGE  PRED  L95  U95; 

/*  LINEARIZING  TRANSFORMATION  MODEL  I  */ 

/*  (TRANSFER  EAGE  TO  l/EAGE)  */ 

DATA  ALL7  ;  SET  ALL1 ; 

IF  EAGE  LE  1100  THEN  DELETE; 

Y= 1 / EAGE ; 

PROC  REG  ;  MODEL  EDTF=  Y; 

OUTPUT  OUT=OUT2  P=  PRED  R=RES  STDI=SEP; 

PROC  PLOT; 

PLOT  EDTF*EAGE ; 

DATA  ALL8  (DROP=  RES  )  ; MERGE  ALL7  OUT2 ;  BY  INDEX; 

PROC  UNIVARIATE; 

DATA  ALL9  (DROP=OBS) ;  SET  ALL 8 ; BY  INDEX; 

W= (Y-0. 000297) * (Y-0 . 000297) ; 

SXX=  0 . 000164*0 . 000164*  (180-1)  ; 

SEPRED  =  112.854*  SQRT (1+ (1/180) + (W/SXX) ) ; 

L95=PRED- 1 . 645*SEPRED; 

U95=PRED+1 . 645*SEPRED; 

PROC  PRINT  DATA=ALL9  ;VAR  INDEX  EAGE  EDTF  PRED  U95  L95  ; 

PROC  PLOT  ; 

PLOT  EDTF*Y= ' A'  PRED*EAGE= ' - '  L95*EAGE= '  .  ' 

U95*EAGE= ' . '  /OVERLAY; 

PLOT  PRED* EAGE = ' * '  L95*EAGE= ' . ' 

U95*EAGE= ' . '  /OVERLAY; 

DATA  EXTEN2 ; 

DO  EAGE=1500 0  TO  35000  BY  100; 

PRED=411. 737+1315827* (l/EAGE) ; 

W= ( (l/EAGE) -0 . 000297) * ( (l/EAGE) -0 . 000297) ; 

SXX=  0.000164*0.000164*  (180-1)  ; 

SEPRED  =  112.854*  SQRT (1+ (1/180) + (W/SXX) ) ; 

L95=PRED-1 . 645*SEPRED; 

U95=PRED+1 . 645*SEPRED; 
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OUTPUT; 

END; 

PROC  PRINT  DATA=EXTEN2 ;  VAR  EAGE  PRED  L95  U95; 

/*  LINEARIZING  TRANSFORMATION  MODEL  II  */ 

/*  (TRANSFER  EDTF  TO  LOG (EDTF) )  */ 

DATA  ALL10;  SET  ALL1 ; 

IF  EAGE  LE  1100  THEN  DELETE; 

EDTFLOG  =  LOG (EDTF); 

PROC  REG  ;  MODEL  EDTFLOG  =  EAGE ; 

OUTPUT  OUT = OUT 3  P=  PRED  R=RES  STDI=SEP; 

PROC  PLOT; 

PLOT  EDTF*EAGE ; 

DATA  ALL11 (DROP=  RES  )  ; MERGE  ALL10  OUT3 ;  BY  INDEX; 

PROC  UNIVARIATE; 

DATA  ALL12 (DROP=OBS) ;  SET  ALL 11 ;BY  INDEX; 

W= (EAGE-4146) * (EAGE-4146) ; 

SXX=  2521553* (180-1)  ; 

SEPRED  =  0.2742  *  SQRT (1+ (1/180) + (W/SXX) ) ; 

L95=PRED- 1 . 645*SEPRED ; 

U95  =  PRED+1 . 64  5*SEPRED ; 

PROC  PRINT  D ATA= ALL 1 2  ;VAR  INDEX  EAGE  EDTF  PRED  L95  U95; 

PROC  PLOT  ; 

PLOT  EDTFLOG * EAGE =' A'  PRED*EAGE= ' - '  L95*EAGE= ' . ' 

U95  *EAGE= '  . '  /OVERLAY; 

PLOT  PRED* EAGE = ' * '  L95*EAGE= ' . ' 

U9  5  * EAGE = '  . '  / OVERLAY ; 

DATA  EXTEN3 ; 

DO  EAGE=14000  TO  20000  BY  100; 

EDTFH=EXP (7 . 3787-0 . 000197*EAGE) ; 

PRED  =7.3787-0. 000197*EAGE; 

W= (EAGE-4146) * (EAGE-4146) ; 

SXX=  2521553* (180-1) ; 

SEPRED  =  0.2742  *  SQRT(1+(1/180)+(W/SXX)); 

L95  =  PRED- 1 . 64  5*SEPRED ; 

U95=PRED+1 . 645*SEPRED ; 

EXPL95=EXP (PRED- 1 . 645*SEPRED) ; 

EXPU95=EXP ( PRED+1 . 645*SEPRED) ; 

OUTPUT ; 

END; 

PROC  PRINT  DATA=EXTEN3 ; VAR  EAGE  PRED  L95  U95  EXPL95 ; 


/*  FLEET  RETURN  DATA  */ 


DATA  ALL 14  ;SET  ALL 2  ; 

IF  DTF  GE  1200  THEN  DELETE; 

PROC  PLOT; 

PLOT  DTF*AGE= ' * ' ; 

PROC  REG  ;  MODEL  DTF=  AGE  ; 

OUTPUT  OUT = OUT 5  P=  PRED  R=RES  STDI=SEP; 
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PROC  PRINT; 

PROC  UNIVARIATE; 

DATA  ALL 15 (DROP=  RES  )  ; MERGE  ALL 14  OUTS;  BY  INDEX; 

PROC  PRINT  DATA  =  ALL15 ;  VAR  INDEX  DTF  AGE  PRED; 
PROC  PRINT  ; 


DATA  ALL16 ;  SET  ALL15 ; BY  INDEX; 
W= (AGE-2821.968)* (AGE-2821.968) ; 


SXX=  1928683* (202-1) ; 

SEPRED  =  61.32  *  SQRT ( 1+ ( 1/2 02 ) + ( W/SXX) ) ; 
L95=PRED- 1 . 64  5* SEPRED ; 

U95  =  PRED+1 . 64  5* SEPRED ; 

AGE3  0  =  3  0 ; 

PROC  PRINT;  VAR  INDEX  AGE  DTF  PRED  U95  L95; 
PROC  PLOT  ; 

PLOT  DTF*AGE='A'  PRED* AGE =  '  -'  L95*AGE=/  .  ' 
U9 5  *AGE= '  .  '  AGE 3 0  * AGE=  7  - ' /OVERLAY ; 

PLOT  PRED* AGE = ' * '  L95*AGE=,./ 

U95*AGE='  .  y  AGE3  0 * AGE=  7  - ' / OVERLAY ; 


DATA  EXTEN5 ; 

DO  AGE  =  10000  TO  16400  BY  100; 

PRED=  823.375-0.0429  *AGE ; 

W= (AGE-2821. 968) * (AGE-2821.968) ; 

SXX=  1928683*  (202-1)  ; 

SEPRED  =  61.32  *  SQRT ( 1+ ( 1/202 )+ (W/SXX) ) ; 
L95  =  PRED- 1 . 64  5  *SEPRED ; 

TEST=PRED- 1 . 699*SEPRED; 

U95=PRED+1 . 645*SEPRED ; 

AGE3  0  =  3  0 ; 


OUTPUT ; 

END; 

PROC  PRINT  DATA=  EXTEN5 ; VAR  AGE  PRED  L95  U95  TEST; 
PROC  PLOT; 

PLOT  L95*AGE='*'  AGE 3 0 *AGE= '  . ' /OVERLAY ; 


/*  WEIGHTED  LEAST  SQUARE  MODEL  */ 

/*  (DTF /AGE  VS  1 / AGE )  */ 

DATA  ALL20  ;SET  ALL2 ; 

IF  DTF  GE  1200  THEN  DELETE; 

NEWDTF=DTF/AGE  ; 

NEWAGE  =  1 / AGE ; 

PROC  REG  ;  MODEL  NEWDTF=  NEWAGE; 

OUTPUT  OUT = OUT 21  P=  PRED  R=RES  STDI=SEP  ; 

PROC  PRINT; 

PROC  UNIVARIATE; 

DATA  ALL 2 2 ; MERGE  ALL2  0  OUT21;BY  INDEX; 

U95  =  PRED+l . 64 5* SEP ; 

L95  =  PRED- 1 . 64  5  * SEP ; 

PROC  PRINT  DATA=ALL2 2 ; VAR  AGE  DTF  PRED  L95  U95  SEP; 
PROC  PLOT; 
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PLOT  RES*  PRE:0='*'  ; 

PROC  PLOT; 

PLOT  L95  *AGE= ' * '  PRED*AGE= ' * '  U9 5 *AGE=' *' /OVERLAY; 

/*  POWER  TRANSFORMATION  MODEL  */ 

/*  (DTF**-3  VS  AGE)  */ 

DATA  ALL23  ; SET  ALL 2 ; 

IF  DTF  GE  1200  THEN  DELETE; 

NEWDTF=DTF*  *  ( -  3 )  ; 

PROC  REG  ;  MODEL  NEWDTF=  AGE; 

OUTPUT  OUT = OUT 2 4  P=  PRED  R=RES  STDI=SEP  ; 

PROC  PRINT; 

PROC  UNIVARIATE; 

DATA  ALL2 5 ; MERGE  ALL23  OUT24;BY  INDEX; 

U95  =  PRED+1 . 645*SEP; 

L95  =  PRED-1 . 645*SEP; 

PROC  PRINT  DATA=ALL2 5 ; VAR  AGE  DTF  PRED  L95  U95  SEP; 

PROC  PLOT; 

PLOT  RES*  PRED= ' * '  ; 

PROC  PLOT; 

PLOT  L95  *AGE= ' * '  PRED* AGE= ' * '  U9 5 *AGE= '*' /OVERLAY ; 


DATA  F1M4EXT; 

DO  AGE=900000  TO  1000000  BY  100; 

PRED=  (1.75*(10**  (-9) ) )  +  (4.87*  (10**(-13) ) )  *AGE ; 

W= (AGE-2821. 968) * (AGE- 2821 .968) ; 

SXX=  1928683*  (202-1)  ; 

SEPRED  =  SQRT(6. 204807* (10** (-19) ) )  *  SQRT (1+ (1/202 )+ (W/SXX) ) ; 

U95= (PRED-1. 645*SEPRED) ** (-1/3) ; 

L95= (PRED+1. 645*SEPRED) ** (-1/3) ; 

AGE3  0  =  3  0 ; 

OUTPUT ; 

END; 

PROC  PRINT  DATA=  F1M4EXT; VAR  AGE  PRED  L95  U95  ; 

PROC  PLOT; 

PLOT  L95*AGE=' *'  PRED*AGE= ' . '  U95*AGE='['  AGE3 0*AGE= '.'/OVERLAY; 


DATA  ONE (  DROP=OBS ) ;  SET  ALL1 ; 

FUME=EAGE/ 365 ;  AGE=EAGE/365 ; 

PROC  PRINT;  VAR  LOTID  AGE  FUME; 

DATA  ONE;  SET  ALL1 ;  IF  AGE  <  3  THEN  DELETE; 


/*  WITHIN- LOT  MODEL  */ 


PROC  MIXED; 

CLASS  LOTID; 

MODEL  FUME  =  LOTID  AGE* LOTID  /  NOINT  S  PREDICTED; 


/*  FIXED  EFFECTS  REGRESSION  MODEL  */ 

PROC  MIXED; 

CLASS  LOTID; 

MODEL  FUME  =  AGE  /  S  PREDICTED; 
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/*  RANDOM  EFFECTS  REGRESSION  MODEL  */ 


PROC  MIXED  ABSOLUTE; 

CLASS  LOTID; 

MODEL  FUME  =  AGE  /  S  PREDICTED; 

RANDOM  INT  AGE  /  CL  G  S  SUB=LOTID  TYPE=UN  ; 
PROC  PRINT; 


/*  COMPUTE  PRED  SEPRED  U95  L95  FOR  RANDOM  EFFECTS  MODEL  */ 


DATA  TWO (DROP=  OBS) ;  SET  ONE;  BY  LOTID; 

PRED=3 . 81021477-0 . 14105394  *  AGE; 

SEPRED=SQRT ( (  0 . 070791) + (2*AGE* (-0 . 00450653) ) 

+  (AGE*AGE*0 . 00030661)  +  (0 . 05926819)  )  ; 
U95=PRED+1 . 645*SEPRED; 

L95=PRED- 1 . 645*SEPRED; 

PREDD=PRED*3  6  5 ; 

U95D=U95*365 ; 

L95D=L95*365 ; 

PROC  PRINT;  VAR  LOTID  AGE  FUME  PRED  PREDD  U95D  L95D; 


/*  GROUP  SHELF-LIFE  ESTIMATION  */ 


DATA  EXTEN7; 

DO  AGE= (5000/365)  TO  (10000/356)  BY  (100/365); 
PRED=3 . 81021477-0 . 14105394  *  AGE; 

SEPRED=SQRT ( (  0 . 070791) + (2*AGE* (-0 . 00450653) ) 
+ ( AGE* AGE* 0 . 00030661) ) ; 
U95=PRED+1 . 645*SEPRED; 

L95=PRED- 1 . 64 5* SEPRED; 

PREDD=PRED*365 ; 

U95D=U95*365 ; 

L95D=L95*365 ; 

AGED=AGE*365 ; 

OUTPUT ; 

END; 

PROC  PRINT;  VAR  AGE  AGED  PRED  PREDD  L95D 


/*  SHRINKAGE  ESTIMATION  */ 


DATA  ' 
CARDS 

THREE;  INPUT  LOTID  INTERCEP  SLOPE 

i 

11142 

-0.06196541 

0.00479259 

11143 

-0.11611040 

0.01007674 

11144 

-0.16769671 

0.01415533 

11145 

-0 . 09516341 

0 . 00494655 

11146 

-0 . 08786149 

0.00599071 

11147 

-0.21248452 

0.01167817 

11148 

-0 . 11397421 

0.00853627 

11149 

-0 . 07768438 

-0.00123725 

11150 

-0 . 15424392 

0.00950736 

11151 

-0 . 12749461 

0.00861269 

11156 

-0 . 15799910 

0 . 00818934 

11157 

-0 . 17588616 

0 . 01418657 

51 


11158  -0.13427054 

11159  -0.23813672 

11160  -0.22564272 

11181  -0.16269825 

11182  -0.14882165 

11183  0.17588730 

11184  0.21035950 

11185  0.44493213 

11188  0.24783741 

11189  0.22284679 

11190  0.18150475 

11191  0.22307183 

11192  0.20898367 

11193  0.19615700 

11194  0.17762001 

11195  0.16493381 

20000  0.00000000 

PROC  SORT;  BY  LOTID; 
PROC  PRINT; 


0.00509910 
0 . 01251881 
0.01415942 
0.00759863 
0.01126227 
-0 . 01074349 
-0.01320753 
-0.02624763 
-0.01470002 
-0.01325530 
-0 . 01253802 
-0.01355682 
-0 . 01170315 
-0 . 01181337 
-0 . 01152574 
-0 . 01078221 
0 . 00000000 


DATA  FOUR;  SET  THREE  ; 

INTERC1=INTERCEP  +  3.81021477; 

SLOPE1=SLOPE-0 . 14105394 ; 

PROC  PRINT; VAR  LOTID  INTERC1  SLOPE1; 

DATA  FIVE;  SET  THREE  ; 

INTERC 1 = INTERCEP  +  3.81021477; 

SLOPEl=SLOPE- 0 . 14105394 ; 

DO  AGE= (5000/365)  TO  (15000/356)  BY  (100/365); 
TSSL  =  INTERC1  +  SLOPE1  *  AGE  ;  BY  LOTID; 
TSSLD=TSSL*365 ; 

AGED=  AGE  *365; 

OUTPUT- 

END; 

PROC  PRINT;  VAR  LOTID  AGE  AGED  TSSL  TSSLD 
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APPENDIX  C  FIGURES 


The  figures  shown  in  this  appendix  are  described  in  Chapter 
II  and  IV  are  as  follows  : 

Figure  C.l  -Figure  C.3  Standard  Linear  Regression  Model  I. 

Figure  C.4  -Figure  C.6  Standard  Linear  Regression  Model  II 

Figure  C.7  -Figure  C.8  Linearizing  Transformation  Model  I. 

Figure  C.9  -Figure  C.ll  Linearizing  Transformation  Model  II 

Figure  C. 12 -Figure  C.14  Standard  Regression  Model  Based  on 

Fleet  return  data. 

Figure  C.15  Weighted  Least  Square  Model. 

Figure  C. 16 -Figure  C.17  Power  Transformation  Model. 

Figure  C. 18 -Figure  C.21  Individual  fitted  Models. 

Figure  C.22  Group  shelf -life  estimation. 
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FUME  TIME  (UNIT  :  DAY) 

400  800  1200 


Figure  C.l.  Fume  time  vs  propellant  age  based  on  the 
standard  linear  regression  model  I. 
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FUME  TIME  (UNIT  :  DAY) 

200  400  600 


Figure  C.2.  Group  shelf-life  estimation  based  on  the 
standard  linear  regression  model  I. 
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Figure  C.3.  Residual  analysis  based 
regression  model  I . 
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on  the  standard  linear 


FUME  TIME  (UNIT  :  DAY) 

400  800  1200 


* 
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Figure  C.4.  Fume  time  vs  propellant  age  based  on  the 
standard  linear  regression  model  II. 
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RESIDUALS 


Figure  C.5.  Residual  analysis  based  on  the  standard  linear 
regression  model  II. 
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Figure  C.6.  Group  shelf -life  estimation  based  on  the 
standard  linear  regression  model  II. 
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Figure  C.7.  Fume  time  vs  propellant  age  based  on  the 
linearizing  transformation  model  1. 
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Figure  C.8.  Residual  analysis  based 
transformation  model  X. 
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Figure  C . 9 .  Fume  time  vs  propellant  age  based  on  the 
linearizing  transformation  model  II. 
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RESIDUALS 


Figure  C.10.  Residual  analysis  based  on  the  linearizing 
transformation  model  II. 
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Figure  C.ll.  Group  shelf-life  estimation  based  on  the 
linearizing  transformation  model  II. 
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Figure  C.12.  Fume  time  vs  propellant  age  based  on  the  fleet 
return  data. 
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Figure  C.13.  Group  shelf-life  estimation  based  on  the  fleet 
return  data 
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Figure  C.14.  Residual  analysis  based  on  the  fleet  return 
data 
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RESIDUALS 


square  model . 
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RESIDUALS 


Figure  C.16.  Residual  analysis  based  on  the  power 
transformation  model. 
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Figure  C.18.  Fitted  within- individual  models  based  on  the 


standard  linear  regression  model  I. 
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Figure  C.19.  Fitted  wi thin- individual  models  based  on  the 
standard  linear  regression  model  II. 
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Figure  C.21.  Fitted  within- individual  models  based  on  the 
linearizing  transformation  model  I. 
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Figure  C.22.  Group  shelf -life  estimation  based  on  the  random 


effects  regression  model. 
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